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Research Paper 25 


Calibration of Meter Line Standards of Length at the 
National Bureau of Standards 


Benjamin L. Page 


The results of intercomparisons of the total lengths of several meter bars and of ¢ali- 


Drations of the 


Information 


which such measurements can be repeated was obtained 
and comments made as to the precautions necessary 


standards of length 


1. Introduction 


It is known that the length of a linear standard, 
at anv specified temperature, cannot be assumed to 
remain the same over a period of years. The changes 
in length may be attributed to secular changes, 
although these changes ure not necessarily uniform 
throughout the length of any linear standard. It 
is essential, therefore, that linear standards be re- 
calibrated at sufficiently frequent intervals so that 
these calibrations may be relied upon to the requisite 
accuracy in their use. 

The standard length in the United States is Proto- 
type Meter 27, kept at the National Bureau of 
Standards. In addition to this meter bar, the Na- 
tional Bureau of Standards has other line standards 
of length that are used as working standards, and 
nine of these were included in the calibrations de- 
scribed in this paper. To calibrate these standards, 
intercomparisons by pairs are made, using different 
combinations in order to determine the most probable 
corrections to their total lengths. As these inter- 
comparisons require measurements of the highest ac- 
curacy, only one pair of bars can be compared on 
any one day. A complete intercomparison, there- 
fore, requires an observing period of months, and 
these complete intercomparisons are made at inter- 
vals of several vears. A previous paper! 
the determination of the total lengths of some of the 
present NBS standards and also the results of that 
and earlier intercomparisons. 

Inasmuch as a complete intercomparison of 10 
meter bars is a time-consuming labor, the question 
arises as to Whether or not statistical methods would 
diminish this labor without significant loss of ac- 
curacy. ‘This paper discusses this question and gives 
the analysis. It also reports the intercomparisons 
that have been made since the report in 1934 and 
shows the trends in the length changes of these bars. 
Redeterminations of the calibration corrections of 
some of the subintervals of the subdivided linear 
standards have also been made. 


describes 


2. Description of meter bars 
Nine meter bars were intercompared in 1952 and 
1953, and the results of these comparisons are dis- 


! Lewis V. Judson and Benjamin L. Page, J 
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subintervals of some of these bars for the past twenty vears are reported. 
both on the degree of stability of 


these standards and on the precision with 
rhe results are critically analyzed 
in the precision comparisons of line 


cussed. Four of these bars are made of platinum- 
iridium, one of 42 percent nickel-steel, three of invar, 
and one of stainless steel (Fe-Cr). 

The four platinum-iridium bars have a modified 
X cross section generally referred to as the Tresca 
section; the others have the more familiar H cross 
section, and all intervals are defined by lines ruled on 
the plane of the neutral axis. When in use, each bar 
is supported at the Airy points, which are located 
28.55 em in each direction from the center of the bar. 

Prototype Meter 27 is the primary standard of 
length of the United States, and meter 21 is a second 
bar of the same series. The other two platinum- 
iridium meters, 4 and 12, are of an earlier melt known 
as Alloy of 1874. The four nickel-steel bars are: 
Meter 39, made of invar, was obtained in 1903; 
meter 153R, made of 42 percent nickel-steel, with an 
inlaid platinum strip on which the graduations are 
ruled, was obtained in 1911; meter 814B, made of 
invar, was obtained in 1931; and meter 752, made of 
a type of invar termed “Fix Invar’’, was obtained in 
1938. The other bar of this group, meter 50, is made 
of stainless and was graduated by C. G. 
Peters, formerly of the Bureau, by use of wave- 
lengths of light. 


steel 


The four nickel-steel bars were obtained from the 
Société Génevoise d’Instruments de Physique of 
Geneva, Switzerland, and are graduated at inter- 
vals of 1 mm throughout their lengths. Meter 4 
is similarly graduated. Meters 27 and 21 have only 
the 1-m interval, which is defined as the interval 
from a central line of a group of three at one end to 
the central line of a similar group at the other end. 
Meter 12 was recently regraduated by C. G. Peters, 
while a member of the Bureau. Originally, this bar 
had rather inferior graduations at the terminal ends. 
This bar is considered as two linear standards for the 
purpose of these calibrations, because two 1-m inter- 
vals were graduated on it with a common zero; one 
interval is nominally correct at 0° C and the other at 
20°C. 

On all the bars subdivided to millimeters there is 
a l-mm interval subdivided to tenths ruled outside 
of the terminal graduations. These may be con- 
veniently used for the determination of the screw 
values of the micrometer microscopes. 

Meter 50 has only a zero line and a 1-m line. 











On each bar two parallel lines approximately 0.2 


mm apart are ruled at right angles to the graduations. 
Only that portion of a graduation lying between the 
two parallel lines is considered when measurements 
are made. These longitudinal lines are also very 
convenient when alining the bars in the comparator 

Table 1 gives coefficients of the linear thermal ex- 
pansion of the meter bars, previously determined 
at the Bureau International des Poids et Mesures 
BIPM) and the National Bureau of Standards 
NBS This table gives the values of a and 6 of 
the general equation 


L.=L,(1 


is the length of the bar at any temperature 
C, and Ly is the length at 0° C 


Ww here L, 


near 20 


| ABLI | (‘oe ficients of erpans on of the mete ha 
a 
Coet! 
Ml 
} All 
‘ Pla iu 1 S.O210 mts 
{ x74 s O14 KX more 
> (1874 1 & HO14 OES 
iT & Holo ~ (Ors 
Inv { 
Soi 4. ckel | ood OOOS 
2 Fix Invar 0. 07 + O40R 
S14B Ir ul 113 
A) S nl 1(Fe-( 10.2 
Meters 27 and 21 are of the principal national series of prototype ete 2 
nd 4 of the allov of 1874° meters 39, L53R, 7 ind 814B have serial numbers of 
the Société d'Instruments {de Physique: and et ) was graduated at the 
National Bureau of Standard Nileter 153R ha ul nl 1 platinur trip 
vhich the graduations are ruled 


3. Description of Longitudinal Comparator 


used in these calibrations 
room in the 
It is mounted 


The comparator (fig. 1 
located in a constant-temperature 
basement of the NBS South Building 
on & massive concrete pier whose upper surface is 
approximately 20 cm below the level of the floor of 
the room and entirely independent of it. This com- 
parator is a longitudinal comparator; the relative 
motion of the microscope beam and of the bars be- 
ing calibrated is in the direction longitudinal to the 
axis of the bars. In the NBS comparator the bars 
are stationary with their axes colinear, and a carriage 


Is 


with a beam bearing two microscopes is movabh 
in a longitudinal direction. The microscopes may) 
be clamped to the beam so that a simultaneous com. 
parison of two bars may be made. Subdivisions of 
a barup toe 2m in length (approximately 80 in.) may 
also be calibrated. The comparator has an over-all 
length of 4' m. 

The longitudinal movement of the 
effected by two electric motors veared to provide 
linear speeds of the carriage from 0.1 to 10 mm/see 
Automatic braking is provided so that the carriage) 
can be stopped within a few microns of any desired 
position. An automatic cutoff is also provided for 
stopping the motion of the carriage to prevent its 
running off the wavs of the bed. The total trave 
of the carriage is 130 em. 

The beam on which the microscopes are mounted is 
attached to the movable carriage through a svsten 
of self-alining ball bearings and a link, so that the! 
beam will have a minimum strain resulting from any 
distortions of the carriage caused by longitudina 
displacement of any magnitude up to 130 cm, re. 
gardless of the distance between microscopes 

The two microscopes have tubes approximately 
15 mm in diameter and 450 mm in length and have 
provision for adjustment. The total magnification 
as used in these calibrations is approximately 16 
diam. The smallest division of the micromete: 
drums is equivalent 0.50 uw, measured the 
focal plane of the microscope. By use of a vernier- 
scale adjacent to the micrometer drum, it is possibk 


carriage js 


to ut 


to estimate to 0.005 u. 

The graduations on the standards are illuminatee 
by means of the vertical illuminators in the micro 
scopes, using low-voltage lamps.  Heat-absorbing 
glass installed between the lamps and microscopes 
greatly reduces the heat radiation from the lamps ane 
therefore the temperature evradients within the bars 
being calibrated. Each microscope is mounted in ¢ 
holder that can be moved along the ways of th 
beam, and may be securely clamped at any position 
Also provided is a slow-motion mechanism for thy 
longitudinal movement of the microscopes. Th 
range in separation of the microscopes is from 1) 
to 210 em. 

Two main supports for the bars are prov ided, am 
each is capable of complete adjustment for focusing 
the bars, so that the parallel longi 


and alining 
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tudinal lines on the bars are parallel to the motion 
of the microscopes and are centrally located in the 
fields of the microscopes Two adjustable roller 
supports, about | em in diameter, are clamped to 
each of the main supports, with the proper sepa- 
ration, for supporting the bars. These supports are 
mounted in the inner of two tanks, which may serve 
as a double bath, the outer one being suitably 
Water may be used in either or both 
tanks and circulated by means of motor-driven 
propellers. Apertures, about 25 mm in diameter, 
are provided in the cover directly above the bar 
graduations that are under observation, 

During the measurements the entire instrument 
is protected against external sources of heat by a 


insulated. 


thick quilted insulating cover, so that only the 
micrometer boxes are exposed, 
4. Calibration of the Meter Bars 
The latest fundamental measurements on this 


group of nine meter bars (considered as 10 for the 


purpose of these calibrations) were begun in 1952 


and completed in 1953. The results of previous 
fundamental and auxiliary measurements of the 
total lengths of these bars, together with some 


repetitions of present measurements of some of the 
subintervals to show the stability of various stand- 
ards over a period of years and the repeatability 
of present comparisons, are included in. this paper. 
Another important part of the work was the cali- 
bration of certain of the subintervals of the bars 
subdivided millimeters to provide accurate 
secondary standards for the calibration of line stand- 
ards that are shorter than 1 m. 


to 


4.1. Determination of Total Length 


All measurements on the bars were made on the 
longitudinal comparator, with the temperature of 
the room thermostatically controlled at 20° ©, The 
bars were allowed to remain in the comparator for 24 
hr, so that they would be in thermal equilibrium at 
a temperature very near 20° C before the measure- 
ments were beeun All measured differences in the 
lengths between any two bars were reduced to 20° C 
by using’ the values for their respective coeflicients 
of linear thermal expansion shown in table 1, before 
computing their respective For 
venience. the corrections to the standards are given 
instead of the actual lengths. The actual length is 
equal to the nominal length plus the correction; 
hence, a plus correction indicates that a standard is 
longer than the nominal leneth. 

Prior the in 1953. the 
thermal conditions of the bars in the comparator 
were Investigated The temperature differences be- 
tween bars, as indicated either by thermocouples or 
by mercury thermometers, were in agreement 
approximately 0.005 deg ©. 

In order to determine the effect that water in the 
outer tank would have on the temperature of the 
bars in the inner tank. the outer tank was filled with 


corrections. con- 


to comparison 1952) and 


to 


water at nearly the same temperature as the bars. 
After the water had been circulated manually, there 
was no perceptible reduction in the temperature 
gradients within each bar nor was the difference in 
temperature between the bars. The motors 
were then used to circulate the water for 2 min. 
Over the period of an hour after the motors had been 
turned off the temperature of the air rose 0.03 deg C, 
and the temperature of the bars remained constant. 
The conclusion is that the thermal conditions in this 
comparator are considerably less stable when the 
circulating bath is used, because of the heat generated 
by the motors. 

The 10 selected for standardization were 
placed in the comparator, two at a time and in a 
random manner with respect to their relative posi- 
tions, either to the right or to the left. To determine 
the magnitude of any possible bias because of the 
location of a standard in either the left or right end 
of the comparator, several auxiliary intercomparisons 
of groups of 5 and 6 of the bars were made. A 
special observing schedule suggested by the Statisti- 
cal Engineering Section of the National Bureau of 
Standards, whereby the bars were mounted in a 
symmetrical manner, showed that there was possibly 
a small bias, which was considerably smaller than the 
reproducibility claimed for the measurements. All 
measurements were made with the bars in air at a 
controlled temperature of approximately 20° C, and 
no water was circulated in the outer tank. 

To obtain the difference in length of two bars, two 
microscopes separated by a distance of 1 m were 
focused on the 0- and l-m graduations of one of the 
Micrometer observations were then made on 
this bar, the procedure being to record 2 observations 
on the left microscope, 4 on the right, and then 2 on 
the left. The carriage was then moved so that the 
microscopes were over the O- and 1-m graduations 
of the other bar, and these graduations were brought 
into focus. Observations were then made on this 
bar. For convenience in computing, the bar at the 
right was always designated as A and the other as B. 
The comparisons between each pair of bars, A and 
4, were repeated 10 times, and the temperature of 
each bar was recorded. The procedure was to ob- 
serve first on B, move the carriage to the right, ob- 
serve on «1, repeat the observations on A, return the 
carriage to the left, and observe on B. This eyele 
was repeated 5 times, and thus a total of 10 differ- 
ences was obtained. The mean of these 10 differ- 
ences was reduced to 20° C, and the difference in the 
length of the bars, L--A-— Bb, thus determined. 
This method was continued until each bar had been 
compared with all others, and a total of [NCNV—1)]/2, 
or 45 differences, had been determined. The results 
obtained are shown in table 2. 

The reduction of these observed differences by 
the method of least squares Was carried out by the 


’ 


use of the form? shown in table 3. 


less 


bars 


bars. 


The differences, +, between the observed and 
calculated differences are shown in table 3. These 
A. Pérard and C. Volet, Les métres prototypes du Bureau International, 


rravaux et Mémoires du Bureau International des Poids et Mesures 21, (1945). 











to 0.20 thre 


The 


0.02 
0.08 


residuals range from nna 


1 


average, neglecting sigh, Is 


error of the caleulated value of the differences in the 
the 


corrections between any two bars, is given by 


formula 


= 


} oe ~ a \ 2 


where >> 
N, the number of observational 


and NV 


equations ; 


the number of unknowns. In the example considered, 





N; is equal to 9 and N equal to 45. For the case 
given in table 3, 7 is calculated to be +0.03. 

r. W. Wright and J. F. Hayford, Th 
137, 143, (1% 


probable 


is the sum of the squares of the residuals; 


During the period from 1933 to 1953, inclusive 
seven series of determinations of the corrections to? 
different groups of the 1G6-m bars under consideration 
These determinations were made with 
vroups of 7. 6. 7 Y . 10 bars All observa- 
tions were made at approximately 20° C, and th 
results redueed 20° C. The reductions of th 
observations were carried out by the method of least 
squares i all The corrections to the total 
lengths of the bars, together with the new corrections 
for the prototype neter determined ut the} 
BIPM and authorized in 1933 by the Eighth General 


were made 
(, 8, ana 


; 


to) 


CUuses 


bars 


International) Conference on Weights and Meas- 
ures, are given in table 4 

The results of any wuxiliary comparisons, mad 
primarily because they served as checks in the in- 
terim between fundamental determinations, are not 


included in the table 

Several points should be noted in considering the 
data in table } The seven values determined lo! 
the corrections to meter 4 are consistently smalle 
than the value reported Dy BIPAI, whereas 6 of thy 


ri values determined for meter 21) and tof the 

values determined for meter 12, before regraduation 
are larger than the values reported by BIPM 
The differs neces between thr eorrections obtained by 


by NBS nnd those reported by BIPNI for prototype 


meters 4, 12. and 21 are +4 0.22, and O.1ld py 
respective ly If we disregard the 1936 and 1988 
Values for meters 12 and 21, respectively, the maxi- 


the NBS corrections is 0.25 wu The 


ud ranee wm 
maximum range in the corrections for meter 4 is 
1) 99 

Figure 2 shows the change in the determined cor- 
rection to the total leneth of each of the two proto 
tvpe mete bars, 4 and 21. from 1933 to 1953, inelu- 
sive, and of meter 12 from 1933 to 1947, inelusive 
and from 1949 to 1953, inclusive 


1 


Krom the above data it appears that there is con 


siderable evidence that these three bars rat have 
changed inh length, two lengthening nnd one shorten- 
ing, by about 0.2) w sinee first received. It is 


rather diffieult to state to what extent the differences 


between the BIPM and NBS values represent 
changes in the actual length of the bars As the 
measurements at the NBS were made under most 
favorable conditions and with extreme care, it Is 
believed that the NBS ts justified in using the latest 
determinations instead of those obtained several 


decades ago, as representing the true value for the 
eorrections to be applied to the respective metel 
length 


bars when calibrating line standards of 
submitted to the Bureau 
Referring again to table 4, it will be noted that 


the five laboratory standards of length, meters 39 
l53R, 752. S14B, and 50, have changed in length 
during the 20-vear period by a much greater amount 
than have the platinum-iridium meters. The thre 
nickel-steel meters, 39, 153R, and S14B, and stain- 
less steel meter 50 have continually elongated during 
this period, but Fix Invar meter 752 has shortened 
The magnitudes of the changes in length between 


successive determinations, however, have generally 
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tended to decrease during the period from 1935 to un 
1953. as shown in figure 3 
, Into 
How long i period might be required for these con 
standards to become sufficiently stable to eliminate few 
frequent standardizations is difficult to estimate the 
These laboratory standards are frequently used \ 
between their restandardization, for measurements by | 
of moderate precision, where the changes in length wit] 
over a comparatively short period of time are not obse 
significant. When they are used in measurements ales 
of the highest precision, they are compared with at squ: 
least one platinum-iridium standard as a check on foot 
their assigned corrections T 
The changes in the corrections during the periods | ones 
between 1933 and 1938, 1938 and 1947, and 1947 | 3 iS a 
—_— - ia S as 
and 1953, are shown in table 5 o \ 
z 
Meter 39 is one of the first laboratory standards ° hela 
‘ : . : w { 
obtained by NBS, and it is of interest to note the | g 
average yearly change in length over a period of | = 
. . rae . Yu 
nearly a half century The changes in length for 
this bar during the period from 1907 to 1953 are 
shown in figure 4. The average changes in length 
per vear, for the first 4, the following 10. and the 
last 32 vears, are 1.2, + 0.6, and 0.2 uw, respec- 
tively. This shows that the magnitude of the 
changes in length is gradually becoming smaller 
The 
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As the work required to observe and determine 15 
for the differences between the lengths of a 
croup of 10 standards is rather laborious and con- 
sumes considerable time, it Was decided to use the 
information obtained in the 1952 and 1953 series of 
comparisons for a study of the possibility of using 


values 


fewer measurements for future intercomparisons of 


the same number of standards. 

According to a scheme? calculations were made 
by using 30, 25, and 15 directly observed differences 
with the remaining differences computed from the 
observed values The reduction of these differences, 


observed and computed, by the method of least 
squares was carried out by the use of forms (see 
footnote 2) shown in tables 6, 7, and 8. 


The scheme for using 30 directly observed differ- 
ences and selecting the remaining 15 to be computed 
is as follows 


Arrange the numbers of the 10-m bars as shown 
below 
27 1 |2a 12b 
27 | 34) 153R 
} 21 72 SI4B 
|\2a 9) 752 50 
12h 153R S14B 50 


Then cancel out the rows and columns that contain 
the number 27; this leaves the numbers 752, 814B, 
and 50 the serial numbers of the bars whose 
differences meter 27 are be computed. 


as 
from 27 to 
Restore the rows and columns and next cancel out 
the rows and columns that contain number 4; this 
then leaves the numbers 39, 153R, and 50 as the 
numbers of the bars whose differences from meter 4 
are to be computed. This procedure is continued 
until each number has been cancelled out of its 
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respective rows and columns and until all the differ- 
ences between the 15 pairs of bars, whose differences 
are to be computed from directly observed differ- 
ences, are determined. ‘To compute these differences, 
the following method, illustrated by 752—27, was 
used. First, the mean of the differences between 
752 and each of the other bars with which it was 
directly compared and for which there is a correspond- 
ing directly observed comparison of that other bar 
with 27 was found. Next, the mean of the differ- 
ences between 27 and each of the others with which 
it was compared, and for which there is a correspond- 
ing directly observed comparison of that other bar 


with 752, was found. Finally, the mean values 
were added algebraically. This is the computed 
value for 752—27, or 179.82 uw, as is shown in 


brackets in table 6. The other bracketed figures 
shown in table 6 were computed in a similar manner. 
All further computations to determine the corrections 
to the bars in this case are the same as in table 3. 

The scheme for selecting and using 25 directly 
observed differences and for leaving the remaining 
20 differences to be computed ts as follows: 


Arrange the numbers of the 10 bars as shown 
below. 

Group | 27 1 l2a 2b 21 

Group 2 39 1538R 752 SHB 50 

The 25 directly observed differences to be used 


are those which are indicated by subtracting each of 
the numbers in group 1 from each of the numbers in 
group 2, i. e., 89—27, 153R—27, 752—27, 814B—27, 
50— 27, 50—21. The remaining 20 differ- 
ences are computed in the same manner that was 
used when 30 directly observed differences were used, 
and are the bracketed figures shown in table 7. All 
further computations to determine the corrections 
to the bars in this case are the same as in table 3. 

The scheme for using 15 directly observed differ- 
ences and computing the remaining 30 differences is 
the opposite of the scheme for which the 30 directly 
observed differences were used and the remaining 15 
computed. The 30 computed differences were ob- 
tained in the same manner that was used where 15 
and 20 differences were determined, and are the 
bracketed figures shown in table 8. All further 
computations to determine the corrections to the 
bars in this case are the same as in table 3. 


In computing the probable errors in tables 6, 7, 
and 8, the same formula (see footnote 3) was used 


asin table 3. This was done for simplicity, although 
the probable errors thus obtained may be somewhat 
too small. Alternative methods for computing the 
differences and probable error are given by Youden 
and Connor (see footnote 4). 

The corrections to the bars obtained in the calibra- 
tion in which all differences were directly observed, 
and those obtained when only 30, 25, and 15 were 
directly observed, are tabulated in table 9. The cor- 
rections to the various bars are in good agreement, 
even though a lesser number of direct differences 


were observed. This will be true only when all the 





observed differences in length between bars are ac- | number (30, 25, and 15) of directly observed differ- 
curately determined ences. The conclusion is that in the comparison of; 

The values in columns 6 to 8 of table 9 were ob- | basic standards, a full calibration, where all the dif- 
tained by taking differences between the correspond- | ferences between bars are directly observed, may be, 
ing values in columns 3 and 2, 4 and 2, and 5 and 2, | justified However, for secondary standards or 
respectively. The average difference in the correc- | others submitted for calibration, a lesser number of 


tions (columns 6 to 8) increases with decrease in the | direct comparisons should be sufficient 
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5. Calibration of the Subintervals 


Methods for the calibration of the subintervals of 
il standard of length have been described nh Bureau 
Circular 329, Calibration of a Divided Seale, by 
Lewis V. Judson. In the recent calibrations of the 
subdivided meter the methods in 
Cireular 329 have been used In certain cases the 
method ° of double or cross calibration has been used. 

The calibration of the subintervals of subdivided 
meter bars 4,39, 153R, 752, and 814B was earried out 
in four steps, namely, the decimeter intervals, the 
centimeter intervals of all of the decimeter intervals, 
the millimeter intervals of the 1- and 92-cm intervals, 
and the 1/10-m intervals at each end of the bars. 

Because of these calibrations, it is possible to com- 
pare lengths from 0.1 mm up to and including 1 m 
with calibrated intervals on any of five 
standards 

The results are expressed as “elements of calibra- 
tion’, that is. the corrections which the subintervals 
will have at a temperature at which the whole interval 
has its nominal length. Thus, it is only necessary to 


bars. deseribed 


these 


distribute) proportionately any newly determined 
total correction to the elements of calibration in 
order to obtain the new corrections to the sub- 


intervals 


5.1. Calibration of the Decimeter Intervals 


The computed values for the several elements of 


calibration of the decimeter intervals of the 5 meter 
bars are given in table 10 
the caleulated values for the elements of calibration 

C. FE. Guillaume, L’ét ivaux et Mémoires du 


illonage des échelles divisées, T1 
Bureau International ’oids et re 


The probable errors of 
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TABLE 10 Elements of calibration of meter bars to decimeters 
153 series of calibrations 
Elements of calibration=corrections to intervals, assuming that the meter bars 
ha heir non il gths. Relative length of interval=nominal length+ek 
mer ealibrat 
Elements of calibration of meter bars 
Int 
i iY 153R 752 S14B 
ba Be “be u A 
| 0. 46 7.36 1. 41 1. 91 0.76 
5) 16 1.39 St) 4.18 ‘4 
0 4). O06 0, 98 51 60 74 
ot { ) 2. 26 409 2. 43 SA 
1 2. 36 1. 61 a) 1. 86 4s 
t 2. 28 2.12 2. 37 1. 87 “4 
0 1. 58 0.06 2 x2 1.78 14 
0 * 4), OS &7 4.84 0. 60 1, 56 
; 1.64 0). 94 2. 82 Ww 1. 25 
0 0 0.0 00 0.00 oo 0.00 
| +0. O2 0.02 0. 02 +0). O2 +0. O2 
Probabl I to to to to to 
| o.0 0.03 0. 0. +0), O45 ou 
Mean of two brations 


are also shown 
formula ' 


They were computed by using the 
> 

} 0.67454 /—» 
VA 


where >)r? is the sum of the squares of the differences 
between the observed and computed values for the 


3) 


elements of calibration. As there are (N—1)(N—2) 
equations of condition, A is [CN —1)(N—2)]/2, the 
value r, is the probable error of any s, and z2 is the 
element of calibration to the first interval, z, is the 


element of calibration to the first two intervals, ., is 
the element of calibration to the first three intervals, 
ete 

That the differences between the observed and 
calculated values for the elements of calibretion are 
relatively small is indicated by the small probable 
errors, and that all of the measurements are appar- 
ently of nearly equal accuracy, regardless of the 
character of the lines on the various bars, is indicated 
by the probable errors being of the same magnitude. 

It is of interest to note how well the values for the 
elements of calibration of the decimeter intervals 
determined during the past 25 vears agree with each 
other, and also to note the differences between the 
mean of these values and those obtained 50 years 
earlier. For this purpose, meter 39 was selected as 
an example, and the values determined in 1903 and 
those determined in 1928, 1931, 1982, 1947, and 1953 
are given in table 11. 

The maximum deviations from the mean of five 
determinations from 1928 to 1953, inclusive, range 
from 0.04 to 0.17 uw, the average being 0.11 x. 

The probable cause of some of the large differences 
between the mean of the 1928 to 1953 determinations 
and those determined in 1903 is that the old Saxton 
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comparator used in 1903 was not entirely adequate | Taste 12. F of calibration of mete 
for the work. It is fortunate, however, that the BIPM and NBS 
values for the elements of calibration thereby deter- ae ae San an a 
mined are in fairly good agreement with those ength. R 
obtained by the use of more recent and better — 
equipment -_ 
5.2. Calibration of the Centimeter Intervals NI mur NBs. 
BIPM Rs tees 
In determining the elements of calibration to the 
centimeter subintervals, the double- or cross-calibra- 
tion method was used, that is, a comparison of each ' ) ) OS 
interval with the others in all possible combinations SO 4 “ rt = s 
namely, the centimeter intervals of the first and sixth to4 ’ . ‘ 
decimeters, those of the second and seventh, the i 
third and eighth, the fourth and ninth, and the fifth to ¢ : | 2s 
and tenth tos ’ i i" 
5.3. Calibration of the Millimeter Intervals : 
The elements of calibration of the millimeter s 
intervals of the centimeter intervals, 0 to 1 em and 
91 to 92 em, were obtained by the same method as ‘ P _ ' 
was used to determine the elements of calibration of to 1. 36 8 
the centimeter subintervals to 10 ‘2 419 rf 
For the sake of brevity the results of the determina- cath ” _ : 
tions of the elements of calibration of the subintervals to | 
of only meter 4 are reported in this paper. The two to 14 0: 98 7 ' 
NBS calibrations for the decimeter, centimeter, and 7 _ ~ - a . ET 
millimeter subintervals on meter 4, together with _ 
those reported in 1907 for the same subintervals by 0 to IS 13 1 29 139 ’ r 0 
BIPM, are shown in table 12 ted > =o ~ a > : 0 
The probable errors for the elements of calibration Ito 21 2.49 2 ae : ah ( 
of any centimeter or millimeter subinterval range 0 to 220 ( - ’ in 4 0 
from 0.05 to 0.07 uw, and were computed by the ~e aap > § “4 : - 9 — 1s : > 
use of the formula / to 254 2.2 1. 2 2.0) 1 s4 22 ‘ ” 
0 to 260 19 H ti { 2 P 
27 j Is i 
r= 0.6745 N= vo um s)he : , 0 
I)yV~I ~ “s 2 | +00 ‘| 4 7 
a ti * | rf | : 
r. W. Wright and J. F. Hayford, 1 sdju 0) Lt 2 (2 st 4 OS v 
ind 142 (1906); Dziobek, Ueber die Ermittelung d rf 0.2 0.27 0.1 
, Ma ibe nach der Methode des Durchschiet \ \ 0 4 ; 24 2.88 2.71 ! 
Ka ! \ \ } n 4,1 (1% ( tu 9 “ " 1 
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limeters at | where =v? is the sum of the squares of the differences 


between the observed and computed values for the 
elements of calibration, NV. the number of subintervals 
into which any considered interval is divided, and r, 


TABI E12 Elements of calibration of meter 4 to mi 


BIPM and NBS—Continued 


I NBS . » > ° ° 
mavi is the probable error of the element of calibration of 
NBS mun NBS ‘ : — 
Inte cane [@nuae| Gare any subinterval, . . 
BIPM ‘BS fron he deviations from the mean of two determina- 
=" tions for all of the elements of calibration for meter 
t range from 0.00 to 0.22 yu, the average being 
. “ 3 ’ 0.06 wu. 
at ‘ 0 La) iM 0.07 0.06 ol 4 — ‘ : — 
5 te 308 ; , o | +.te 02 . "3 lhe discrepancies between the values originally re- 
ote ah : 5 - 6 - ms ported by BIPM for these elements of calibration 
and the values determined by NBS vary from —0.3 
; { Jl 21 21 in) 4 2 = . a . . 
0 te _— 0 { Is L 2] 03 1 9 to +O.) w. However, over a period of years repeated 
Oto 480 3 r : 1B 03 4 determinations of the elements of calibration of the 
1H) +. 2 +. 24 +. 36 2 +. 2 decimeter subintervals on meter 4 have shown highly 
70 ’ 2 “4 9 0 consistent results, and it is indicated, therefore, that 
ost 3 .. 3 ~ 2 . v8 the NBS values are very nearly the present true 
1 to 500 2 2s 2. 36 Os values for all of the elements of calibration for the 
( | as ' 1 33 0 2 . ont a ° ° . 
subintervals. ‘Taking into consideration that the 
20 ‘ - . oa = s NBS values were determined by using a comparator 
10 0 2.71 2.87 272 05 +3 of more modern design, with more accurate auxiliary 
~~ > La . 7 02 equipment and under more favorable conditions 
aa , :' than were available in 1907 at BIPM. it is believed 
si) 2.1 2.8 2s 2.82 r2 that the NBS is justified in adopting these new values 
«00 ; 22 > 35 2 2s 07 2 as representing more nearly the true values than 
610 ' | ; | dais those originally determined. 
630 24 oe 2. 22 2. 14 5.4. Calibration of the One-tenth Millimeter 
pel : = ; Intervals 
toa | 4 es = : The ‘o-mm subintervals at each end of the meter 
ye 3 2 8 > 0 > 82 02 bars were also calibrated by using the same double- 
‘ : ) #8 0.76 1 ‘tg or cross-calibration method used in the calibration 
of the centimeter and millimeter subintervals. The 
j " 9 9 } 
" { { 0.0 0 i 
4 i at ye L i t. 1 TABLE 13 Elements of calibration of the one-tenth millimete 
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NBS 1953 series of calibrations 
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values determined for the elements of calibration for 
meter 4 are shown in table 15 

The probable error for the elements of calibration 
The devia- 


of any ‘io-mm subinterval is + 0.04 u 
tions from the mean of two sets of determinations 
range from 0.00 to 0.10 uy, the average deviation 


being 0.03 u 
6. Conclusions 

In making the measurements on the meter bars, 
the results of which are reported in this paper, every 
reasonable precaution was taken so that the results 
would reflect the accuracy anal precision that is 
obtainable with the present NBS equipment 

In the comparison of the bars to determine their 
total corrections and the corrections to the subinter- 
vals, the smallness of the residuals indicates that the 
standards were in good thermal equilibrium while 
the observations being made. Although the 
observations were made when the temperatures of 
the bars were very close to 20° C, the small residuals 
also indicate that the values used for their respective 
coefficients of linear expansion are very nearly the 
true values. In calibrating these meter bars, the 
effect of any error of 0.01 deg C in temperature varies 
from approximately 0.01 u/m for the invar bars to 
0.10 u/m for the stainless steel (Fe-Cr) bat 

To a large extent the effect of the imperfections 
of the lines and surfaces of the bars was compensated 
for by sufficient repetitions of the observations. — It 
is believed that the final correction to the total length 
of each of the standards calibrated is not in error by 
more than 0.2 uw, and in most cases is not in error by 


were 


more than 0.1 uz. 
There is no evidence that the lengths of the plati- 
num-iridium bars have changed a significant amount 


14 


during the past 20 vears, whereas there is definite 
evidence that the other bars have changed. These 
changes, however, have become smaller during the 
ensuing Vears 

Regarding the possibility of making fewer measure- 
ments when intercomparing a group of standards, 
the 1952 and 1953 series of calibrations showed re- 
markably good agreement with the final corrections 
to the bars if a lesser number of direct differences had 
been considered. This could only be the Cuse if the 
observed differences inh length between the standards 
are all in excellent agreement In the ealibration of 
basic standards, a full calibration should made 
where all differences between lengths of the bars are 
However, where standards of a 
calibrated, a 


directly observed. 
lower degree of accuracy are 
number of direet comparisons would suffice 

As it is known that bars of ans materials are 
unstable in length, the 1952 and 1953 series of com- 
parisons continues the program of maintaining all 
line standards of length and of obtaining additional 
knowledge about them, so as to enable the NBS to 
meet the needs and demands of the present and future 


lesser 


for increased accuracy, 


The author expresses his appreciation to Lewis 
V. Judson, Chief, NBS Length Section, for his untir- 
ing assistance and guidance and for the many helpful 
suggestions through the vears during which the data 
reported in this paper were obtained, and to John S. 
Beers for his assistance with reeent observations and 
computations 


WaAsHINGTON, July 30, 1954. 


exa 
mail 
not 
fror 
reg 
doe 


Seri 


hite 
ese 


the 


ure- 
rds, 
re- 
ions 
had 
the 
irds 
n of 
ade 
are 
f a 


sser 


are 
OMl- 
all 
nal 
5 to 


Lure 


‘WIS 
itir- 
pful 
lata 
ns. 


and 


yf Research Paper 25¢ 


Enthalpy and Specific Heat of Four Corrosion-Resistant 
Alloys at High Temperatures 
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tvpe 446, containing 26 percent of chromium 
67 percent of nickel and 30 percent of copper 


Pour 


measured with a drop method and a precision Bunsen ice calorimeter, 


fo 
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1. Introduction 


During the past several vears the specific heats of 
a large number of substances have been measured 
precisely between 0° and 900° C at the National 
Bureau of Standards In each case this involved 
measurements not only with the sample in a suitable 
container, but also with the empty container as well, 
in order to determine by difference the net effect of 
the sample alone. Each such container consisted of 
a thin-walled evlinder, in. in diameter, and was 
usually composed of one of four alloys: 80 Ni-20 Cr, 
Monel, or stainless steel type 347 or 446 (Am. Lron & 
Steel Inst. designation). These alloys have numerous 
uses because of their relatively high resistance to 
corrosion, but apparently few measurements of their 
specific heats at elevated temperatures have been 
reported. 

In the experiments recorded in this paper, the 
specimens used were themselves empty containers 
of the type previously used for measurements on 
other materials, such aluminum oxide. The 
accuracy is somewhat less than in the earlier work 
because there is no cancellation of the small heat 
losses resulting from subtraction of blank measure- 
ments with an empty container. In addition, the 
following two factors should be pointed out, which, in 
the present case, somewhat limit the accuracy with 
which such measurements on only a few specimens 
can be applied to all specimens of the same general 
allov: (a) As these four alloys are relatively hard 
materials, the previous thermal and mechanical 
history of a given specimen may sometimes cause it 
to exist in some less stable state that 
somewhat different specific heat; (b) a single name is 
often applied to a group of alloys whose chemical 
compositions, although nearly the same, are not 
exactly so. Fortunately, the specific heats of the 
main constituent elements of these four alloys are 
not very different, thus tending to reduce any error 
from this source. In faet, only in the temperature 
region of a transition or where two phases coexist 
dloes it seem likely that these factors may lead to 
serious lack of reproducibility in specific heat. 
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Discontinuities in specific heat were found 


and the stainless steel tvpe 446 


2. Experimental Procedure 


Very detailed descriptions and evaluations of the 
method and apparatus used for measuring the 
enthalpy are contained in two papers [1] ' and in a 
later paper [2], describing the modified form of 
apparatus used in the present work. In brief, the 
method consisted of heating the specimen in a 
furnace to each of several known temperatures and 
then dropping it into an ice calorimeter, by means 
of which the heat evolved in cooling the specimen 
to 0° C was measured. The specimen container 
(which was the actual specimen in the present 
work) was suspended in the furnace and calorimeter 
by a No. 32 wire of 80 Ni-20 Cr. This container 
and wire were connected by a short platinum- 
rhodium wire to which were attached two thin 
platinum disks that also entered the calorimeter. 
These disks greatly reduced upward loss of heat 
from the specimen. The helium inside the container, 
as well as a gold gasket in the case of 80 Ni-20 Cr, 
also entered the calorimeter during the drop. Helium 
was continuously passed through the cores of the 
calorimeter and the furnace, in order to reduce 
oxidation of the sample and to hasten the attain- 
ment of thermal equilibrium. The temperature of 
the furnace was kept constant during each run to 
within +0.01 deg for a sufficient length of time to 
allow the sample to attain the same temperature. 
The temperatures were measured up to 600° C 
with a platinum resistance thermometer and above 
600° C with a platinum-—platinum-rhodium thermo- 
couple, both of which had been calibrated in terms 
of the International Temperature Scale of 1948. 
The calibration factor of the ice calorimeter, which 
was based on numerous earlier electrical determina- 
tions at the Bureau, was taken to be 64.646 calories 
per gram of mercury (1 cal=4.1840 absolute joules). 

Because, as mentioned above, these alloys were 
not themselves in containers, it was necessary to 
estimate and apply two small but appreciable 
corrections to each heat measurement, namely, for 


Figures in brackets indicate the literature references at the end of this paper. 
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a) the heat lost during the drop through the region 
between the furnace and the calorimeter, and (b) 
the heat delivered by the irrelevant parts of the 
system dropping into the calorimeter. The value 
used for the first correction was the mean of the caleu- 
lated and observed values previously discussed [2] 

The four alloys used were of commercial grade 
In those cases where measurements were made 
with several specimens of the same alloy, the speci- 
mens were taken from adjacent parts of the same 
sample of rod or tubing. Except in the case of the 
Monel. which was stated to be in cold-drawn con- 
dition, they were heated to above 1,000° C and al- 
lowed to cool in air before beginning measurements 
of their enthalpy. Kach specimen had 2 mass 
between 8 and 17 g, and carried with it into the calo- 
rimeter 0.3 to 0.5 g of platinum metals. The S80 
Ni-20 Cr was accompanied also by 0.3 0.6 2 
of gold, and the other three alloys, by about 0.1 ¢ of 
80 Ni-20 Cr. In the course of the measurements, 
all the alloys except the Monel acquired oxide 
coatings in varying amounts up to a few hundredths 
The thermal contributions of the small 
irrelevant materials, which 
stitute correction (b) of the preceding paragraph, 
were estimated from = specific-heat data in the 
literature and were corrected for by deduction from 
the results found. 

A separate specimen of each alloy, taken from the 
same sample as those used in the thermal measure- 
The resulting 


to 


of a gram. 


amounts of these con- 


ments, was analyzed chemically 
compositions are given in table 1, and fall within 
the ranges of nominal! cOMPposi tion for alloys of 


these designations [3]. 


3. Results 


The enthalpy measurements were on two specl- 
mens of 8O Ni-20 Cr, on one each of Monel and 
stainless steel 446, and on six of stainless steel 347 
of which two were not run above 450° C Several 
duplicate were made on each 
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wre 4 
men at each furnace temperature, the uverage 
deviation from the mean being 0.08 eal @>'. or 
about 0.05 percent of the total heat measured The 
different specimens of stainless steel 347 showed 
small differences that were systematic with tem- 
perature and averaged 0.25 percent of the measured 
heat The corrected net heat values for the two 
specimens of St) Ni-20 Cr. however. differed hy 


approximately 1.5 percent at each furnace tempera- 
ture. As this difference is many times the precision 
of the measurements, it is attributed to differences 
in the two specimens 

Relative to that at O° C, the mean 
enthalpy at each furnace temperature, corrected as 
deseribed above but unsmoothed. Is viven for each 


observed 


alloy in table 2 The mean values for the two 
specimens of SO) Ni-20 Cr are listed separately 
because of their considerable differences. In the 
last two columns are given (in calories per gram of 
allov) the major corrections that were made to the 
values for stainless steel 446. The corresponding 


corrections for the other three alloys were very 
similar in magnitude, ; 

The following equations, derived from the values 
of table 2, represent the relationship between 
enthalpy (in eal ¢g and temperature (in deg C 
The average departure of the observed values from 


these equations Is about 0.02 eal ¢ less than 0.1 
percent The equations for SO Ni-20 Cr were 
derived from the mean values observed with the 
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H,-—Il O.1033¢ sS(107-°)t (O° to 150° © (1 
H.-W 0.33 + 0.10798 
2 32(10~>°)t (150° to 600° C (2 
HT,—ll 1.37 +-0.11471 
2.31(10°°)4 600° to 900° C (3 
Vonel: 
H—ll 0.1010f+2.20(107°)f 0° to 300° C). (4 
Stainle ss steel 347 “4 
H,—I1,=0.12508t+ 1.703(10~°)¢ 
9.31 low 1.) (¢#+ 273.2)/273.2 0° to 900° C 5 
Stainless steel 446: 
Hi—H 0.10781#+ 5.371107) 0° to 500° C). (6 
He—H 6.83 4+-O.16178¢ 
23 O000(10 600° to 900° C). (7) 


Differentiation of these equations leads to the fol- 


lowing equations for the specific heat (in cal g 
dee C 
SO Ni-20 Cr 

Hy 0.1033 + 7.60107) 0” te 150° © S 

C,—0.1079+ 4.64(107°)f 150° to S00° ©). q 

c 0.1147+4.62(10°°)f HOO? to 900° C 10 
Vonel: 

Hy 0.1010+ 4.40610 f 0° to 300° CC 1] 
Stainless steel 347: 

( O.12514+3.41010°°)¢ 

$.04/(t4-273.2 0° to 900° C 2 

Stainless steel 446: 

$s; 0.1078+-1.074(107*)# 0° to 500° C 13 

e 0.1618 4-2209(10 600° to 900°C 14 


Specific-heat values calculated for even temperatures 
from eq (Ss 14) are listed in table 3. 

The specific heats of the four allovs are shown 
graphically in figures 1 to 4. The lines represent 
the smooth values vlven by eq (Ss to 14 and the 
points correspond to the unsmoothed corrected values 
calculated for short temperature intervals from pairs 
of successive values given in table 2 

The representation of the specific heat of 80 Ni-20 
Cr by straight lines in figure 1 is convenient and 
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TABLE 3 Specific heat (smoothed) 
Specific heat, cal g-' deg C 
remperature Stainless steel 
80 Ni-20 C Monel 
347 446 
( 
0 0. 1033 0. 1009 0.1104 0. 1078 
2 1052 1021 1123 1105 
“« 1071 1033 1142 1132 
100 1109 1054 1176 1185 
a) 1171 1097 1233 1203 
$00) 1217 1142 1283 1401 
100) 1264 1326 1508 
un) 1310 1370 163 
600 14 1410 21 
700 1470 1448 1695 
SOO S15 . 1487 1631 
900 1563 1525 1620 


within the accuracy of the data, which however do 
not justify the assumption of a real discontinuity 
in slope at 150° C. On the other hand, real dis- 
continuities do appear for 80 Ni-20 Cr near 600° C 
(fig. 1) and for stainless steel 446 between 500° and 
approximately 650° C (fig. 4). Although the total 
enthalpy changes over these regions were measured 
as precisely as in the other cases, the exact variation 
of specific heat with temperature in these regions of 
transition was not established because the enthalpy 
measurements did not cover sufficiently small tem- 
perature intervals. Otherwise the smoothed specific 
heats represented by eq (8 to 14), table 3, and the 
lines in figures 1 to 4 are believed to be accurate 
within +2 percent, at least for the actual specimens 
measured. 
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4. Discussion 


The following is a comparison of the specific-heat 
values inh cal ty deg % obtained in the present 
work with those given in reference [3]. A value of 
0.107 given for SO Ni-20 Cr [3, p. 1060], presumably 
near room temperature, is the same as that given for 
50° Cin table 3. For stainless steel 347, the averag 
between 0° and LOO0° C is given as 0.120 [3, p. 566] 
the present work gives a corresponding value 5 per- 
cent lowe! For wrought stainless steel 446, 0.144 
at 100° C is listed [3, p. 555]; the corresponding value 
of this paper for stainless steel 446 is about 20 percent 
lower. A value of 0.127 at 1,260° © is listed for 
standard Monel [3, p. 1049], but is not necessarily In 
disagreement with the present results, which extend 
up to only 300° C. Compared with a value for Monel 
of 0.0997 at 20° CC. obtained with an adiabatic eal- 
orimeter [4], the present result is 2 percent higher 

A comparison of the specific heat found for stain- 
less steel 347 with corresponding values measured al 
the National Physical Laboratory (London) [5] fora 
stainless steel of somewhat similar composition (8.14 
percent of nickel, 19.11 percent of chromium, 0.37 
percent of manganese, 0.68 percent of silicon, 0.08 
percent of carbon, and 0.60 percent of tungsten) may 
be significant. The specific heats of the latter steel 
recorded in table 4 are the average values found for 
successive 50-deg intervals whose mean temperatures 
were those listed. The values given for stainless steel 
347 were calculated from eq (12). Though the regu- 
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absence of a transition in the case of stainless steel 
347, the specific heat of the other steel shows, between 
550° and 600° CC, a small anomalous that 
ordinarily attributed to some sort of transition. 

The present investigation shows, over a tempera- 
',a very marked increase in 
the specific heat of the particular sample of stainless 
steel 446 investigated. Though the details of this 
transition may depend somewhat critically on the 
exact chemical composition of the sample, it appears 
that the effect observed is to be associated with the 
magnetic-transformation point atl approximately 
600° C that has been reported for iron-chromium 
allovs ol this composition 16] 

Several have found an anomalous 
behavior in a number of physical properties of 80 
Ni-20 Cr in the temperature region just below and 
above 600 * An abnormal change with tempera- 
ture of the thermal expansion from 500° to 600° C 
was found [7], and a marked increase in the specific 
heat as the temperature rises above 500 
reported by others [8]. Recent measurements have 
been made of the specific electrical resistivity and 
specific heat of this alloy over the temperature range 
from = 400 760 C’ {9}. These measurements 
show, with Increasing temperature, a decrease of 
electrical resistivity and an anomalous increase in 
specific heat, tentatively explained by assuming 
that the alloy contained regions of very-short-range 
order whose number Or size decreases iis the tem per- 
ature The data indicated furthermore 
that the specifie heat was somewhat less after the 
alloy had been quenched from 1,000° C than after it 
had been heated at 500° C and then cooled slow ly. 

The discontinuity in the specific heat of SO Ni-20 
Cr near 600° C shown in figure 1 represents a behavior 
not encountered in most transitions. As the temper- 
ature the specific heat rather abruptly shifts 
to a higher, parallel line 
the measurements, all the heat absorbed by the alloy 
Is represented by the area under the curves Since 
a “hump” such as that often found in the specific- 
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heat curve of a substance in its transition region is 
missing, there is no heat of transition. In addition, 
more recent unpublished results obtained at the 
Bureau, which are considerably more precise than 
those on 80 Ni-20 Cr, unmistakably show similar 
behavior. These recent values are for two alloys 
manufactured under the trade name Inconel (approx- 
imately 8 percent of iron, 76 percent of nickel, 15 
percent of chromium, and 0.05 percent of carbon), 
in which the amounts of the two major constituents 
are not greatly different from the amounts in 80 
Ni-20 Cr. 

The behavior of these alloys may be analogous to 
that of cobalt. On the basis of enthalpy measure- 
ments by a number of investigators, Kelley [10] gave 
the heat of its transition near 400° C as zero. This 
seems plausible, because according to X-ray analysis 
the transition is from a close-packed-hexagonal to a 
face-centered-cubic phase [11], a process that involves 
no change in atomic coordination number and only a 
slight change in the interatomic distance. It has 
similarly been found that the nickel-chromium alloys 
with less than 45 percent of chromium have a face- 
centered-cubie structure (presumably at room tem- 
perature) [12], and an explanation similar to that 
suggested for cobalt may hold in the case of 80 
Ni-20 Cr. 


Most of the experimental work was supported by 
the Atomic Energy Commission. D. C. Ginnings, 
Anne F. Ball, and Ann W. Harman, of the National 
Bureau of Standards, assisted in making a large 
number of the measurements. The chemical anal- 
vses of the 80 Ni-20 Cr, Monel, and stainless steel 
347 were performed by B. B. Bendigo, also of NBS. 
Stainless steel 446 and its chemical analysis were 
supplied by the Allegheny Ludlum Steel Corporation. 
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Effect of Temperature on the Tensile Properties of a 
Commercial and a High-Purity 
70-Percent-Nickel-30-Percent-Copper Alloy 


William D. Jenkins, Thomas G. Digges, and Carl R. Johnson 


Shor ne tensile tests were made at temperatures ranging from 75° to 1,700° F on 
wo 70-percent-nickel—30-percent-copper alloys The experimental evidence showed that 
lisconti < flow occurred in specitnens of both alloys fractured in tension at temperatures 
ranging ft 00° to 1700 1 This phenomenon was attributed to strain-aging at the 
lower te eratures and to reerystallization accompanied by grain growth at the higher 
Cerri pe i Variations in chemical composition of the two alloys also affected the degree 

rain-nei strength and ductility properties and fracture characteristics 


1. Introduction 


Tensile tests were made 
from 75 1,700° F 
nickel »O-percent-copper allovs of the same initial 


erain size and structure but differing in purity. One 


at temperatures 


to on two lots of 70-percent- 


ranging 


alloy was of relatively high purity, whereas the other 


was a commercial grade that contained appreciable 
amounts of manganese, tron, and cobalt. The tensile 
properties of the high-purity alloy as evaluated at a 
lower range of temperatures 320 212° F) have 
been investigated 1] and a study of its creep he- 
how Ih progress 


to 
havior is 
2. Previous Investigations 


The influence of the mechanical 
properties of commercial grades of 70-percent-nickel 


temperature on 


0-percent-copper alloys (Monel) has been discussed 
n several publications |2, 3], but little information 
s available on the effeet of impurities on these 


properties. An inspection of the results given in the 
tables |2] shows that the tensile strength of hot-rolled 
Monel increased as the test temperature was raised 
to 600° F without a corresponding decrease 
n elongation However, in the range of 1,000 
1.300° F there is a definite increase in brittleness, 
accompanied by a loss of strength. MeVetty and 
Mochel [4] reported a precipitous drop in values for 
tensile strength and duetility of hot-rolled Monel as 
the test temperature from 750° to 
Q30° KF 

Some creep properties of Monel have been reported 
by Bennett and MeAdam [5], MeAdam, Geil, and 
Woodard [6], Grant and Bueklin [7], Shahinian and 
Lane {8}, and Betty, Eiselstein, and Huston [9]. In- 
dications of the onset of intererystalline fracture, 
aging, and discontinuity in’ temperature-ductility 
relations are evident in some of the results presented. 


from 500 


to 


Wiis increased 


I ure ty Ket | te tl terature reference the end of 
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3. Alloys and Testing Procedures 


The chemical composition of the alloys investi- 
vated is given in table 1. 

The charges were melted in an induction furnace, 
and each allov was cast into a 14- by 14- by 60-in. 
ingot. The ingots were milled, then forged to 8- 
by S8-in. blooms, hot-rolled to 24-in. square billets, 
and then hot-rolled to %-in. rounds. The round 
bars were annealed at 1,100° F. for 16 hr, then at 
1,000° F for 8 hr. The alloys, as supplied by the 
manufacturer in the above condition, had an average 
erain diameter of 0.025 mm. The structures of these 
allovs are shown in figure 19, A and B. 

Tensile specimens, 0.505-in. diameter and with 
2.0-in. reduced sections, were machined from one 
bar of each alloy. Each specimen was heated in 
air in an automatically controlled electric furnace 
and held at the desired temperature for approxi- 
mately 1 hr before testing in tension; the specimen 


Was maintained within +3 deg F of the desired 
temperature. The specimens were tested in a 
hvdraulic-tvpe machine, with the specimen and 
TABLE 1 Chemical composition (percentage by we ight of the 
alloys used 
Chemical composition 
DD 
N Cu Mn F« Si ( 
Hlig! yur 0. OS YY 71 0.01 0.01 0.12 0.017 
( 6. OT 29.79 SY 1. 5 7 15 
S Ca Co O: Ne He 
H purity 0. 002 ( 0. 001 0. OO15 0. 0003 
Commercial 006 0.017 0. 38 OOO! None 0001 


Faint trace 


Not detected 





two trolled to produce a rate of strain of about 1 percent 
minute 

Hardness tests and metallographic examinations 
were carried out at room temperature on the frac. | 
tured) specimens in with 
described procedures [10] 


adapter occupying about 38 ft between self- 
alining holders. The movement of the head during | a 
the test was recorded by means of a “strain 
recorde! of the Templin type Except for SOTLLC 
specimens tested at 600° F, as noted in tables 2 and 
a. the movement of the head of the machine was con- 


stress 


accordance previous| 
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4. Results and Discussion 
4.1. Mechanical Properties 


The experimental results are summarized in tables 
2 and 3 and in figures 1 to 22. A comparison of the 
flow characteristics in tension of the alloys at room 
temperature is obtained by observing the relative 
positions of the true-stress-true-strain * graphs in 
figure 1. The flow stresses at equal strains, stresses 
at maximum load (1/7) and at beginning of fracture 
R,), and rates of work-hardening were greater for 
the commercial than for the high-purity alloy. 
However, the true strains at maximum load and at 
the beginning of fracture and fracture strength * were 
correspondingly greater for the high-purity alloy. 

Photographs of many of the autographic stress- 
strain records obtained on specimens of the alloys 
tested at different temperatures are reproduced in 
figures 2, 3, and 4. The rate of work-hardening 
varied with the amount of strain, the test tempera- 
ture, and composition of the alloys. At the lower 
test’ temperatures, the slopes of the curves were 
relatively steep from the onset of plastic deformation 
to the maximum load. Precipitous drops in the 
curves in the region of maximum load to fracture are 
associated with necking of the specimens, whereas 
the flatness of some of the curves in this region is an 
indication of uniformity in deformation. At equal 
temperatures and strains, the rate of work hard- 
ening of the commercial alloy was usually somewhat 
higher than that of the high-purity alloy. Discon- 
tinuous flow occurred in many of the specimens as is 
indicated by the serrations in these curves. These 
serrations were ofien obtained both before and afte 
the attainment of the maximum load on specimens 
of each alloy. Although the trend was for the serra- 


rt 1 ‘ was determined | lividing the current load by the current 

nimum cross-sectional area of the specimetr rhe true strain was determined 

t sl logarithm o {tl nal cross-sectional area of the speci 
i o the curre m 1u ctional area i 

Frectu ength was determined by dividing the final load by the minimun 
ect ila of the fractured sp en at room temperature 


tions to attain a maximum in the specimens tested 
at 600° F, they were also observed in specimens 
tested at temperatures ranging from 300° to 900° 
F, and were considerably more noticeable in the 
stress-strain curves obtained on specimens of the 
high-purity alloy (fig. 2) than in the curves for the 
commercial alloy (fig. 3). Furthermore, -increasing 
the speed of testing (table 2) by a factor of about 15 
did not materially decrease the magnitude of serra- 
tions in the curves for the specimens of the high- 
purity alloy tested at 600° F (compare curves A and 
(’, fig. 2). The serrations are believed to be asso- 
ciated with strain-aging and other atomic rearrange- 
ments. 

Discontinuous flow was obtained in specimens of 
both alloys tested at the higher temperatures, as is 
illustrated in figure 4. The breaks in the stress- 
strain curves in the region of strains corresponding 
approximately from the drop of beam to nearly com- 
plete fracture are interpreted as being associated with 
recrystallization * and grain growth (fig. 20, E and 
F; fig. 21, C and D; fig. 22, E and F). 

The effect of temperature on the yield strengths 
of the two alloys is shown in figures 5, 6, and 7. 
For the specimens tested at the usual rate, there 
appeared to be a continuous decrease in yield strength 
with increase in test temperature, with the exception 
of the specimens of high-purity alloy tested within 
the range of about 600° to 700° F (fig. 5) The 
slopes of these curves were continuously changing, 
however, and the magnitude of the change was 
affected significantly by the composition of the 
materials; vield strengths of each alloy at 1,500° or 
1.700° F were approximately equal. Increasing the 
speed of application of load (points designated B) 
at 600° F raised the vield strengths slightly. 

Certain investigators have indicated that a linear 
relationship exists between the vield strength and 
the reciprocal of the absolute (Kelvin) temperature 
when the experimental results are plotted on a linear 
or semilog (log of temperature) basis. The results 
obtained in the present investigation are plotted in 
these manners in figures 6 and 7, respectively. 
These assumptions were found to be approximately 
valid for limited ranges of temperature, but for each 
alloy there is a considerable difference between the 
slopes of the curves at high and at low temperatures. 
Moreover, at equal temperatures, there is also a 
difference in the slopes of the curves in each of the 
diagrams for the two alloys. This suggests that the 
activation energy for flow is not a simple function 
of the applied stress, but rather a complex function 
involving stress, temperature, and the size of the 
flow units activated by temperature. 

The tensile strength of each alloy decreased with 
an increase in test temperature from 75° to 300° F 
(fig. 8). However, reversals were obtained in the 


Metals Handbook, 1948, p. 1049), the commercial alloy 
md 1,110° F for 10- and )-percent reduction of area, 
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tensile strength-temperature curves between 300 
and 500° F for the commercial alloy, and between 
500° and 600° F for the high-purity alloy. The 


tensile strength of each alloy decreased with a further 
Between 1.000° to 


increase in test temperature, 

1,700° F, the tensile strengths of both alloys were 
nearly the same. The reversals are attributed to 
strain-aging, and the superior strength of the 


commercial alloy at the lower temperatures is con- 
sidered as being due to impurities restricting the 
mean free path of the flow units within the grains 
Increasing the speed of testing at 600° F tended 
to decrease the tensile strength of the high-purity 
alloy and to increase slightly the tensile strength of 
the commercial alloy. These observations are also 
consistent with the that the impurities, 
acting as barriers, tended to increase the resistance 
to flow of the commercial alloy at the higher speeds. 
The effect of temperature on the fracture strength 
(see footnote 3) is shown in figure 9. It should be 
pointed out that these values were affected by the 
necking characteristic (rim effect of the alloys; 
the rim effect was more pronounced in the high- 
purity alloy than inthe commercial alloy. The frac- 
ture strength of the high-purity alloy attained a 
maximum at 212° F. At the usual speed of testing, 
the strength then decreased as the test temperature 


concept 





700° F. 
Increase 


and 
in 


a change in temperature from 600° to 
again decreased rapidly with further 
temperature. The maximum 
the commercial alloy occurred at a higher temper- 
ature range (300° to 500° F 
trend was about the same as for the high-purity alloy 
At any test temperature below 1,100 j. 
the fracture strength of the high-purity alloy was 
higher than that of the commercial alloy Th 
fracture strength-temperature curves for the two 
allovs nearly coincide at temperatures ranging fron 
1,200° to 1,700° F. 


selected 


fracture strength of? 


_and then the veneral | 


} 


Increasing the speed of testing from the norma 


to the intermediate rate used at 600° F poimts B 


resulted in raising the fracture streneth of the com- 


mercial allov, whereas the strength of the high-purity 
alloy was not materially affected by this change 
However, the strength of the latter allov increased 
markedly as the speed of testing was further in- 
creased compare points Band C 


The relations between true stress and true strain! 


at maximum load and at fracture are shown in figures 
10 and 11, respectively. Strain-aging effects ar 
again evident in the lower range of test temperatures 
Although the trend was for the true strains at mavi- 
fracture the true 
and the were in- 


to decrease as 


temperatures 


mum load and at 


stresses decreased, 
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alloy and from 700° to about 900° F for the com- 
mercial alloy, both curves for the latter alloy show 
reversals at a temperature of about 900° F. 

The values obtained for strain at maximum load 
and elongation and reduction of area (tables 2 and 3 
were not significantly affected by variation in test 
temperature from about 75° to 700° F (fig. 12); the 
values for strain and elongation for each alloy were 
somewhat higher at 75° F than the corresponding 
values at 212° F. The minima obtained in the 
ductility-temperature curves for specimens of the 
commercial alloy at a temperature of about 900° F 
are associated with the brittle behavior of this alloy 
at this temperature. 

Variations in the range of speed of testing used at 
600° F had no material effect on the ductility of the 
allovs. 

The results of hardness tests made at room temper- 
ature on the specimens fractured at different temper- 
atures are summarized in figures 13 and 14. The 
hardness increased continuously with the amount of 
deformation for the specimens of the high-purity 
alloy fractured at temperatures up to and including 


1 000° F and also for the specimens of the commercial 
alloy at temperatures ranging from 75° to 800° F 
(fig. 13 With a further increase in temperature to 
about 1,200° F, the hardness at room temperature 
of each alloy was not significantly affected by the 
amount of plastic deformation. At test, tempera- 
tures above 1,200° F, the hardness generally de- 
creased as the amount of deformation increased; no 
marked change in hardness occurred in the specimen 
of the high-purity alloy, tested at 1,500° F with 
change in reduction in area from about 20 to 42 per- 
cent. The hardness—-reduction-in-area curves for the 
specimens tested at 600° F are at a somewhat higher 
level than those for specimens tested at lower or 
higher temperatures. The curves of figure 14 also 
show that the maximum hardness of each alloy was 
obtained in the specimens fractured at 600° F, 
These curves (figs. 13 and 14) verify the fact that 
combination of strain-hardening and strain-aging 
predominated in the tensile carried out at 
temperatures ranging from about 75° to 800° or 
1,000° F, whereas recovery and_recrystallization 
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predominated at temperatures in excess of 1,200° F, 
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42. Fracture Characteristics and Microstructures 


The effect of tem tut 


erature on the tendenes of the 


loys to contract locally (neck) ts illustrated by the 
hotographs of the fractured specimens (fie. 15 
j At the lower range of temperature, specimens of 
both allovs contracted loeally and fractured in a 
retile manner with a relatively large rim effect 
fig. 16 At the higher range of temperature, the 


specimens deformed more uniformly over the entire 
reduced section and fractured in a relatively brittle 


manner without a pronounced rim effect (figs. 15 
and 16 However, the trend toward brittle be- 
havior at the higher test temperatures was con- 


siderably more prominent for the commercial than 
lor the high-purity alloy. This feature is also shown 
va comparison of the curves of figure 17. It 
noteworthy that the reduction-in-area——distance- 
lrom-fracture curve for the specimen of high-purity 
alloy tested at 1,700° F all the curves for 
specimens of this alloy tested at lower temperatures. 
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fension at different te m peratures, 
il alloy at the right of the photograph 
I i y I J, 1,200° I k, 1,350 " LL. 1,500 F: M, 1,700 k 


This fact again illustrates the complexity of the 
mechanisms of flow and fracture of the alloys at 
elevated temperatures. 

The appearance of the surface and the types of 
fractures obtained in the allovs tested at various tem- 
peratures are shown in figure 18. No surface cracking 
was observed in any of the specimens of the high- 
purity alloy tested at temperatures up to and in- 
cluding 900° F (fig. 18,D); shear-type fractures and 
circumferential strain markings were characteristic 
of these specimens. Except for fewer strain mark- 
ings, these characteristics were also obtained in 
specimens of the commercial alloy fractured at tem- 
peratures up to and including 800° F (fig. 18,A and 
B). Surface were observed in all of the 
specimens fractured at the higher temperatures, and 
the brittleness increased with temperatures (fig. 
I8,D and FE, C and F); the surface appearance was 
similar for the specimens of both alloys fractured 
at 1,200° F (not shown). Thus, the tendency to- 
ward brittle behavior was accentuated by the in- 


cracks 








creased amount of impurities and by an increase im 
test temperatures. 

The influences of temperature on the microstruc- 
tures and on internal cracking of the specimens 
are illustrated in figures 19 to 22, inclusive 

The initially equiaxed grains of the annealed 
alloys (fig. 19,A and B) were elongated in the direec- 
tion of the applied stresses in the specimens fractured 
at temperatures below the recrystallization range 
(fig. 19); the extent of the elongation apparently de- 
creased as the temperature was increased. The 
microcracks were confined to the region of complete 
fracture in the specimens of the high-purity alloy 
tested at temperatures up to and including 900 Kr 
(fig. 19,C and EK) and in the specimens of the com- 
mercial alloy tested at temperatures up to 600° F 
(fig. 19,D). Appreciable microcracking, extending 
from the axis to the surface, occurred in the specimen 
of the commercial alloy fractured at 900° F (fig. 


19.F and H). The formation of an oxide film (top 
edge of the photomicrograph of fig. 19,H) may have 
also contributed to the brittle behavior of this 


specimen 

Numerous microcracks were observed both in the 
region of fracture and at some distance from it in 
all the specimens fractured at temperatures in ex- 





cess of the recrystallization temperatures (fies, 29 
21, and 22). The number of these cracks and th 
extent of cracking, however, varied with the tem. 
perature, composition of the allov, and distance frop? 
the region of complete fracture. At 1,200° F, th 
two allovs were nearly alike with regard to crackin 
(fig. 20.A and B; fig. 22,A and B), and no detectable 
erain growth occurred, f 
At 1,500° F, the microcracks observed in the regio 
of fracture were possibly somewhat less numeroy 
but of larger size in the specimen of the high-purity} 
alloy than in the specimen of the commercial allo 
fig. 20.C and D: fig. 21,A and B This was try 
also in the region about 0.5 in. from the fracture (fig 
22C and D The microcracking wes confined! 


principally to the grain boundaries, and no abnorm; 


grain growth was detected in’ the reervstalliz 
structure } 

At 1,700° F, intererystalline cracking was promi-| 
nent in the vicinity of fracture and as some distany 
from fracture in the specimens of both alloys (fig 
20,E and F; fig. 21,C and D; fig. 22,E and F Strai 
markings 2nd twins were present in the recrystallize 
structures at the time of failure. Appreciable grai 
growth also occurred at this temperature. 
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5. Summary 


Tensile tests were made at temperatures ranging 
from 75° to 1,700° F on specimens of both high-purity 
and commercial 70-percent-nickel-30-percent-copper 
alloys, each with the same initial grain size and 
structure. 

Serrated 
specimens tested at temperatures ranging from 3500 
900° F. These serrations were more noticeable 
for the high-purity alloy than for the commercial 
alloy. Discontinuous flow occurred in 
mens of both alloys tested at temperatures above 
900° F. The serrations were interpreted as being 
associated with strain-aging and other atomic rear- 
rangements, and the discontinuous flow at the higher 
temperatures as being associated with reerystalliza- 
Strain-aging was also mani- 


stress-strain curves were obtained on 


to 


also spect- 


tion and grain growth. 
fested in the strength-temperature relationships of 
both alloys; it was especially prominent in the curves 
for tensile and fracture strengths. Variations in 
composition of the two allovs markedly affected its 
strength properties (yield, tensile, and fracture) and 
reduction of area at temperatures below 1,100° F. 
Above 1,200° F, the strength properties for the two 
allovs were quite similar, but the values obtained for 
ductility (reduction of area, strain at maximum load, 
and elongation) were influenced by the composition 
of the alloys; at 1,200° F, the ductilities of both alloys 
were nearly the same. 

The rates of work-hardening varied with the test 
temperature, amount of strain, and composition of 
the alloys 

The results of hardness tests made at 
perature on specimens of both alloys fractured in 
tension at different temperatures indicated the pre- 
dominance of strain-aging and work-hardening at 
the low range in temperature, a balance of work- 
hardening and recovery in the intermediate range, 
and the predominance of recrystallization and re- 
covery at the higher temperatures used in the tensile 
The hardness at room temperature of the 


room tem- 


; 
ests 


fractured specimens was also dependent 


composition of the alloy and 
deformation 
At temperatures below 900 


both alloys contracted locally 


a ductile manner with a rim 
goo? F, 

uniformly 
manner without a 
and fracture characteristics we 


pronounces 


composition of the alloys 

The final structures and the 
in the specimens fractured at 
varied with the test temperat 
deformation, and compos'tion 


upon the 
the amount of plasti 


F, the specimens of 
before fracturing jy 


At temperatures aboy 


specimens of both allovs deformed mor 
before fracturing in a relatively 


britth 
“it rim. The necking 
re also affected by the 


microcracks produce 


different temperatures 


ure, amount of plast 
of the alloys. 
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The System Lime-Water at 21° C and High Pressures 
Charles E. Weir 
Compression measurements on the binary system Ca(OQH).-H.O at 21° C and high 
pressures show a large discontinuity in volume at 2,800 atmospheres which involves both 
; components l.vidence is presented that the transition is the result of formation of a 
hydrate of lime, Ca(OH)o-n H,O, where 4<n<6 \ pressure composition diagram for the 
system Ca(OH)o-H.O at 21° © is given. 
} 
1. Introduction to obtain a uniform mixture and then stirred with 
Varsol prior to compression. The dry Ca(OH), 
In a recent study of the properties of portland | Was treated in a stmilar manner except that no water 
cements subjected to high hydrostatic pressure, a | Was added. Following compression of this speci- 
large volume discontinuity was noted in systems | Men, the volume was determined by hydrostatic 
containing cement and free water 11] By elimina- weighing in Varsol. The resulting calculated density 
tion procedures it was finally concluded that the | for the Ca(OH). was used to caleulate the volumes of 
volume discontinuity arose from the svstem | Ca(OH), in all experiments in which water-Ca(OH). 
Ca(OH).-H.O. The “present report deseribes the | mixtures were used. The volume of the water was 
results of direct studies on the svstem Ca(OH).-H.O | caleulated from known data for the density of 
at 21° C and high pressures . oo water. The amount of Ca(OH). dissolved was 
' The literature on hvdrate formation in alkaline neglected in this calculation. 
earth hvdroxides is extensive and will not be re- 
viewed here. It is sufficient to note that although 3. Results and Discussion 
numerous reports have been made of hydrates of 
Ca(OH). at room temperature, none appear to be The results of compression measurements on the 
accepted as reliable [2] and no hydrate of Ca(OH), is | various mixtures are given in table 1. The mole 
believed to exist. As a result of the present studies | fraction of Ca(OH), in each mixture is tabulated 
fit is concluded that a hydrate of Ca(OH), is stable | at the head of each corresponding column of com- 
at pressures above 3,000 atm at 21° C, and a rough | pression data. The compressions are expressed in 
pressure-composition diagram is proposed. terms of —AV/\y, where —AV is reckoned from 
the arbitrary zero at 2,000 atm, and Vo is the volume 
2. Experimental Method of the specimen at 1 atm. Values of compression 
at pressures below 2,000 atm, therefore, appear as 


Techniques and apparatus used here for studies of 
compression at high pressures have been deseribed in 
detail [3]. Briefly, an experiment consists in foreing 
a leakproof piston into the bore of a heav -walled 
pressure vessel that contains the specimen immersed 
ina light petroleum distillate (Varsol).  Measure- 
ments are made of the depth of penetration of the 
piston at each 1,000 atm, the pressure being de- 
termined by means of a calibrated manganin pressure 
immersed in the Varsol inside the pressure 

(nother experiment performed with a steel 
bar of comparable volume replacing the specimen 
permits calculation the compression of the 
specimen in terms of that of steel, which is known. 

The specimens used in these studies were mixtures 


rage 


vesse] 


of 


| of powdered Ca(OH , ana distilled water contained 


in a evlindrical steel tube open to the confining 
liquid at one end Ca(OH). was of chemical reagent 


crade dnd was weighed out in the tube. Freshly 
distilled water was also weighed into the tube. so 
weights of both components were known. 
Rach water-Ca(OH), mixture was stirred thoroughly 
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The lowest pressure attained in these 
studies was 1,143 atm rather than 1,000 atm, and 
an additional measurement was made at 2,470 atm 
in most experiments, as shown in table 1. At the 
higher water contents it was not possible to attain 
the maximum pressure of 10,000 atm because of the 
excessive piston travel arising from the large transi- 
tions. In these experiments the initial readings 
were made at 9,000 atm. In experiments with 0.12 
and 0.022 mole fractions of Ca(OH), attempts were 
made to attain 10,000 atm, but the immediate 
crystallization of ice-VI, initiated at a pressure of 
about 9,300 atm, produced too large a piston travel. 
In the measurement with 0.080 mole fraction of 
Ca(OH), the maximum pressure was only 9,100 atm 
and no ice-VI crystallized, but the measurement at 
9,000 atm undoubtedly represents data for a me- 
tastable system. 

The rather complex behavior of the tabular 
data is simplified by considering the compression- 
pressure curves shown in figures 1 and 2. Figure 1 
contains the data for dry Ca(OH), and the mixtures 
of high Ca(QH), content, and corresponding data 
for water and the mixtures of low Ca(OH), content 
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are given in fivure 2 The data lol Hhqguia watel 
were reported by Adams [4] 
The following observations may be made from 
the curves of figures | and 2 
1). The large discontinuity in volume in the 
svstem CalOH -H.O involves both components 
because ho comparable effeet occurs in water or 
drv Ca(OH), alone 
2). Quantitatively, the amount of the diseon- 


tinuity reaches a maximum for some intermediate 
composition of the mixture. 
) evidence exists for formation 
at low water contents, but at high 
ice-V I separates at approximately S.S00 atm 1 
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gmount that increases with increasing water content 
ne polati 


} It mav be coneluded, therefore that 
liquid Watel ‘is present ni pressures above » ON in the 
atm inmixturesricher in Ca(OH), thanapproximately) (SUM 
20 mole pereent This conelusion follows in part HO 
from observations (3) above and in part from th As 
very low compressibilities slopes of the AV/V yp Hon, 
versus P curves) noted above 3,000 atm, whi 
are characteristic of solid phases 

5 In the presence of excess water, th (tiscon 
tinuitv in volume at 2.800 atm is relatively sharp, Where 
but exeess Ca(OH appears to “spread out’ the) lieved 
discontinuity over a wider pressure rane Soli HO : 
Ca(OH therefore appears to act as an impurity comp) 


and produces rounding of the corners of the dis- | water of hydration merely for convenience as it is 


continuity As the phases present above 3,000 atm recognized that this material may be a form con- 
are believed to be all solids in mixtures rich in | taining a complex aquo-cation [6]. As evidence 
y Ca(OH), contamination may be interpreted in | against the former type of compound, it might be 
terms of solid solutions. anticipated that a material containing ordinary 


6). The compressibilities of the products formed | water of erystallization would be less dense and 
above 2,800 atm in mixtures rich in Ca(OH), seem | probably more compressible than the parent anhy- 
to be slightly less than the compressibility of the | drous material. The possibility that the volume dis- 


—— 


parent dry Ca(OH), (see table 2). It must be continuity arises from transformation of amorphous 
noted, however, that in this region two components | Ca(OH), to the more dense crystalline Ca(OH), is 
are most probably present discounted because liquid water would not disappear 


~_ 


7). The freezing pressure of water in mixtures | in such a process and reversibility would not be 
rich in water Is definitely below 9 OOO atm. Inas- expected, 
much as pure water freezes under these conditions The value of » in eq (1) may be bounded in the 
at approximately 8.800 atm [5], the material formed | following manner: In these mixtures the quantity of 
by the transition at 2,800 atm can be only sparingly | product formed will be proportional to the quantity 
soluble of the component present in deficiency. That is, in 

8). The presence of excess water is clearly shown | the presence of excess Ca(OH), the amount of reac- 
by the curvature noted above 3,000 atm in the | tion 1 will be determined by the quantity of water 
curves of figure 2 and by the increasing resemblance | present, whereas in the presence of excess water the 
of these curves to that for pure water as the water | Ca(OH), becomes the determining factor. If the 
‘content increases. The strong dissimilarity be- | discontinuity in volume is taken to be a measure of 
tween these curves above 3,000 atm and those of | the extent of reaction 1, a plot of mole fraction versus 


— 


~~ 


figure | is to be noted volume change at discontinuity should produce two 
el 9). The sharp drop in compressibility noticed | lines intersecting at the composition corresponding 
on formation of the high-pressure phase in both | to the pure compound. Such a graph of the data 
nat figures | and 2 may be taken to indieate disappear- | of table 2 is shown in figure 3. In figure 3, it is noted 
ance of water coincident with formation of the high- | that in the region of excess H,O good linearity is ob- 
| pressure phase tained but in the region of excess Ca(OH), curvature 
appears to exist. This may arise from interaction 
Tawny 2. Volume discontinuity and compressibility data between the components in the presence of excess 
Ca(OH)., as mentioned in observation (5). As a 
aud result of this curvature and the large slope of the 
Mole fraction in volume» Compressibility linear section, it is not possible to place narrow 
— wot fede phase limits on n. From figure 3 it would appear that it 
can be stated with assurance that 3<n<7 and with 

considerable less certainty that 4<n<6. 
From the foregoing considerations, a rough phase 
diagram for the system Ca(OH),-H,O at 21° can be 
fF constructed. This diagram is given in figure 4, and 
- is based on the previous discussion. In the solution 
} 23 ‘ region, however, no measurements have been made, 


as it was believed that no conclusions could be 














The compressibilities and discontinuities of volume Pc : ; . , , 
. ri > 
referred to above are tabulated in table 2 These 7 
1 . a 
lata were all obtained graphically and for this reason 0" 
are not to be considered as precise In addition, the S of ‘ 
} ° ° ° ~ 
f evaluation of the discontinuity in volume is subject a 
, . . : : ¢ 2 0.10 
tent, 0 obvious extrapolation errors incurred in extra- @ “a ; 
. - ¢ 
+ no! Polating to 2.800 atm Che large curvature existing “ Ps 
> ° at 
» yg) in the presence of excess water prevented graphical ~ ; 
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ately ¢CSUumation of compressibilities for mixtures rich in 5 / 
+ HO 3 a 
part ) > ya 
i the) Asa result of the foregoing observations the reac- Deed JA 7 
V/V tion, e "ii 
> y, 
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scorl 6 ai 
os ; : , 0 l l i l l l l 
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10 a 7 It must be emphasized that the data given here 
e S ee = | were obtained on decreasing the pressure, and figur 
t represents these findings. On increasing pressure 
8} the low-pressure forms may persist at considerably } 
higher pressures than those indicated in figure 4 fo 
. the transition lines. This feature is quite commo, 
3 6 in most high-pressure transitions \ 
2 The results reported here may be of significance / 
¥ S reologically because hvdrolvzable calcium salts In 
- minerals subjected to pressures above 3,000 atm may 
& tend to exist in the form of the high-pressure hydrat ) 
3 This transition was originally observed in_ wet 
y hvdrated cement paste and in wet hvdrated tri- 7 
é calcium. silicate, both of which contain Ca(OH), 
However, it does not occur in gvpsum-water mixtures} 
: at pressures as high as 10,000 atm In view of th 
. . a (Ct 5 0.6 —_ 8 ? O anomalous solubility behav lor of calcium compounds Re 
WATER and the large decrease in volume at the transition j an 
Fiaure 4 ee . msrecsure-comneniises gram fo sade however, there appears a distinet possibility ot the the o 
. Ca(OH)»-HLO at 21° ¢ occurrence of such processes at still higher pressures} ‘wed. 
in the earth shia 
formed because of the slight solubility of Ca(OH), naril 
In this region. therefore. the course of the curve has 4. References senen 
been estimated. The initially increasing solubility heen 
tt i " ' 77s _— ll C. bk. Weir and C. M. Hunt (in preparation). 
of Ca(OH), with increasing pressure is inferred from ae it @. Madea: Pics ehmeleall ateceeeas ek Wi et 
the fact that its solubility increases with decreasing pounds, vol. 1, p. 248 (Clarendon Press, Oxford, 1950).{ energ 
temperature 7| Either this solubility decreases as 31 C. EF. Weir, J. Research NBS 45, 468 (1950) RP2160 ) tratic 
is indicated at higher pressures or dissolved Ca(OH 1 LL. H. Adams, J Am Chem soc 53, 3709 1931 hy es 
; 51 International Critical Tables, vol. I\ \icGraw-H 
is transformed into the compound of eq (1), which ts Book Co.. Inc. New York. X. ¥.. 1928 I ean | 
sparingly soluble, since the saturated solution freezes | i¢) J] J. Emeleus and J. 8S. Anderson, Modern aspects meth 
at approximately the same pressure as pure HO inorganie chemistry, p. 191 (D. Van Nostrand Co.) from 
In the latter case, the transition line at 2,800 atm Inc., New York, N. Y., 1952 I ; sion | 
would extend across the solution region Above a6 Lae seo and L. 8. Wells, J. Research NBS 62, 75 (1954 Re 
3,000 atm, the solid phase may be a mixture involving re lin wl 
a solid solution. In the light of the present knowl- mixe 
edge, satisfactors representation of any solid solution brom 
field is of dubious value. Noattempt to include such disk 
a field on the diagram was made WASHINGTON, September 15, 1954 in di 
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Absorption of Radiant Energy by Solid Particles 
in Suspension 


James E. 


\ simple 
SUS} nded ina tr: 
sUspel ded particles 
\pparent deviations 
experimental work 


Increases from Beer’s 


difficulties i 


1. Introduction 


Radiation-absorption measurements have long 
been used for quantitative analysis of materials in 
the gas and liquid phases. Solid materials, however, 
have been troublesome even for qualitative determi- 
nations because of the inherent difficulties in pre- 
paring specimens in a reproducible form. Until 
recently the standard method for study of solids has 
been to suspend finely divided particles in a medium 
such as mineral oil in order to reduce the amount of 
energy lost by scattering. In principle, the concen- 


> tration of such a mixture, or mull, can be controlled 


by careful weighing, and the thickness of the mixture 
ean be kept fixed by proper cell design [1] But the 
method inconvenient, and its usefulness suffers 
from interference by absorption bands of the suspen- 
sion medium 

Recently a new technique [2, 3] has been proposed 
in which finely ground particles of a solid sample are 
mixed with an alkali halide powder, usually potassium 
bromide, placed Ina die, and pressed until a sintered 
disk is formed. At the Bureau, small disks, 9.5 mm 
in diameter and approximately 100 mg in weight, 
are produced hh a simple lie patterned after one 
deseribed by Anderson [4]. Pressure is applied with 
an ordinary shop vise. No special precautions are 
taken for drving the potassium bromice or evacuating 
the system while pressing. Satisfactory qualitative 
infrared spectra have obtained with these 
pellets, and the method seems promising for quanti- 
tative Measurements as well 

The purpose of this report is to present a simple 
theoretical derivation of an absorption law for these 
disks and to infer from the law certain precautions 
to be taken when preparing samples for quantitative 
The discussion applies equally well 
to powders suspended in liquid media.  Mathe- 
matical descriptions of radiation seattered from 
spherical particles have long been known. But they 
are quite complicated and cannot be approximated 
satisfactorily when the particles are about equal in 
size the wavelength of incident radiation. In 
addition, index of refraction curves in the neighbor- 
hood and interior of infrared absorption bands are 


Is 


been 


measurements 


to 


not known for many covalent compounds. <A few 

published dispersion curves 15, 6] show that over 

many Infrared absorption bands refractive index 
Figure t i nad t th tu erence it the end his p 


4l 


theorv 1s developed to describe 
insparent medium and also by an absorbing medium containing transparent 
It is found that in both cases transmittance increases as particle size 


are discus c 1. 


Stewart 


the absorption of radiant energy by particles 


law and other effects that might introduce 


changes very little. Price and Tetlow [7] have made 
the same observation from their work on Christiansen 
windows in infrared spectra. Accordingly, the pres- 
ent study is concerned with the case of no interface 
refraction, requiring that a perfect index of refraction 
match between particle and surrounding medium be 
assumed. In reality this condition can be only 
approximately true at best. 


2. Derivation of an Absorption Law 


Consider first a single particle of material with 
absorption coefficient @ defined by Beer’s law, 


I/T, 


exp(—azr) for a uniform layer of thickness z 
(see fig. 1 do, is an increment of the area s; 
on which the incident plane wave falls with ir- 


radiance ? J), n is a unit vector normal to the particle 
surface, j is a unit vector perpendicular to the radia- 
tion wave front such that the wave travels in the —) 
direction, and d/ is the energy leaving the particle 





through do. Beer’s law can be written for a 
evlindrical increment of volume 
dE =I yn-jdoy exp(yo— mya. (1) 
Radiant flux | ul urface are 
+ 
- r 
i \ 
i i \ 
pono 
\ fo = = === === arias ? > 
ihe fe apap tet do ( 
- } x Y2,Z22 Ki, Vi 2y) - ae Y 
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I irfticle in 


pi opagate ad in the 


Figure 1 adiant-e ne ray Si ld, the ene rgy be ing 


) direction. 








The projected cross-sectional area of the particle on 


the plane of the wave front is 


I] 


or spherical particles of radius 7 the integrations in 


2) are easily done to give the equivalent expressions 
Ti] | Jar lLjexp Par)|/2a*) 
exp 2a 2a llexp Yar) /2a 
2a) 
& ya » 9 
Z Th Xp Zu) 
y Vf y 4 
hr,ayexp 2a) » 
Jones [8S] has investigated the effeet produced in 


the apparent absorption coefficient by holes in a 
laver of absorbing substance. The next step in the 
present development follows his method. Con- 
sider a disk containing masses m, and of sample 
and supporting medium of densities and py. 
respectively The sample particles are assumed to 
be spheres of radius r, and be distributed unt- 
formly throughout the unpressed disk on a cubie 
lattice of constant A. The distance between layers 
of particles parallel to the wave front will be reduced 
by the pressing. It will be evident from what 
follows that the distribution of particles within 
lavers can be random so long as the number of 
particles per layer is constant. In fact, randomness 
within layers will improve the validity of an approxi- 
mation which will be made. The volume of disk 
occupied by sample and by the suspending substance 

' p, and vr py, respectively The number 
of particles in the disk is 


Mi 


pp 


to 


Mi Mi 


Is / 


and the mean distance between particle eenters is 


\ } / iL | Consider first one lave of particles 
in cells of side X. The ratio of the cross-sectional 
area, o,, Of the sample particle contained in a single 
cubie unit to the cross-sectional area, of the 
lattice unit is 
a T) yy rin p 
| 
0 - } i, p ii, p 
and the transmittance for one layer is given by 
I] | : t : fir.da) exp 2a 7 
( o Oo | ' 
Now the number of lavers in the disk is .\ S, 


where S is the disk cross-sectional area. Substituting 


42 


from above, front, 
, | om Mn Mi : j Nf : 

A *( —+— } 6) | =eXP 

rs itp, p p Kq 

radius 


Equation (5) gives the average transmittance of one j gnall 
laver of particles. We assume that the J given by (5 
for one laver is the /,) for the next laver. The irra. 


diance after the radiation has traversed V lavers js 


therefore 
| In 
2a) |: 7 
( 


Il  ( l ; )+" h(r,@aexp 
| 


After the work reported here was completed, a| 
similar discussion appeared in the literature {9 
contamimng an expression for the transmittance of 


) the e 


(discal 


Wem 


potassium bromide pellet in the form J/Jo=« n) | series 
in Which F is the average transmittance of a single | *"™4 
particle and @ contains mass, density, and dimension |?" 
parameters. A comparison of this expression ex-| 
panded in a Maclaurin’s series with eq (7) expande 
in a binomial series demonstrates the equivalence of 
the equations, at least for small particles r| 
Now let us turn to a similar problem in which th wo 
particles are transparent and the suspension medium | YP 
absorbs energy. A situation such as this is encoun- 8 
tered in conventional mulling techniques of sample | _e 
preparation but will not arise for pellets unless Tiewer 
becomes necessary tO Use suspension powders that peer 
have absorption bands. Let @ be the absorption — 
coefficient of the suspension medium end all other| rs ; 
parameters remain eas defined previously. Proceed-|- “y 
ing as above, the average transmittance for a singk pond 
lattice unit is that | 
-_ The ¢ 
J] l | |. \ } idle Ny wre f pared 
a | Sie ) | transi 
; small 
, )( | ad )+ Tp =ra—I 1) g) | taimin 
) 6 o 2a?) j must | 
It h 


The transmittance for the first laver of cells equals} tion b 
the same value and, making the same assumption! |s redt 


used earlier regarding the average J, for the next} have | 
laver, the transmittance of the specimen ts an Un 
ever, 
Th gy} Will re 
sorptic 

3. Discussion a 

band 
It is not clear how eq (2) might be generalize = 
further in order to demonstrate the effect of particle for th. 
shape.” However, some easily integrable shapes cal) jp; ’ 
be investigated. For example, a disk-shaped particle)» wit 
with volume and cross-sectional area equal to thal spectr 
of a sphere with radius r must have thickness 47/3. pein, 
With its flat surface parallel to the radiation wav the oan 


| ferent 


I one 


\ 
irra- 


Crs js 


od 


y 


| of a 


single 
Hsion 


front, the particle would have a transmittance of 
[To )aise = @XP | far/3) For the sphere (J/Jo) spnere 
| =exp( for sufficiently small ar. 
Equation (7) for spherical particles involves the 
radius in a complicated manner. If ar is sufficiently 
small so that terms beyond those of first degree in 


2ar ) 


| the expansions of exp(2ar) and exp(—2ar) can be 
{ disearded, the expression becomes 

o ra] Da 
| Ini / V In | (- ) (14 = ya Par) 

Oo oO ») 
{ , 0 ta) 

\ In ( |- ): 

| oO m 


We may now expand this expression in a Maclaurin’s 
| series, keeping only the leading term, substitute (4 
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fand (6), and make the simplifving assumption 

}m,<m, (Which is true in practice). The result is 
lim In J// am ,/Sp,. 

[his is recognized as Beer’s law for a mass m, of 


fsample in the form of a homogeneous condensed 
film. For larger particles the sensitivity of (7) to 
changes in r can best be demonstrated by assuming 
jsome typical values for the other parameters and 
plotting ///) versus r 2). There is indeed a 
significant dependence on particle SIZe, which is 
most serious for bands with highest absorption co- 
efficient. As these are the bands most likely to be 
ised for quantitative purposes, the necessity for 
controlling particle size is clear. It is also apparent 
that band shape should be affected by particle SIZe. 
The curves of figure 2 predict that a sample pre- 
pared with large particles will have a higher apparent 
transmittance than a sample with the same mass of 
small particles. This applies also to specimens con- 
taining clusters of small particles. In this case 7 
must be taken as the cluster radius. 

It has been observed experimentally that absorp- 
tion bands do become more intense as particle size 
sreduced. Some observations on polyvinyl chloride 
have been made. It is recognized that scatter plays 
an undetermined role in these measurements. How- 
ever, it is not likely that reduction of particle size 
will result consistently in greater scatter at an ab- 
sorption band, and hence in greater apparent absorp- 
tion. Furthermore, no asymmetric variations in 
band shapes, such as Christiansen windows, were 
observed. 

Calculations from eq (7), using typical values 
lor the parameters, provided the encouraging result 
[that the apparent transmittance is independent of 
m, Within the practical measuring limits of modern 
spectrophotometers. In order to verify this’ ex- 
perimentally, two pellets were prepared containing 
the same quantity (2.6 mg) of stearic acid but dif- 
ferent amounts of potassium bromide (42.25 mg 
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using the 
0.714em 2:4 


calculated from eq 


0.224 v2: p 1 @: pi 0.2 g/em NS 


and 294.85 mg Transmittance was measured rel- 
ative to base lines drawn through maxima adjacent 
to the minimum in question. The transmittances 
of corresponding bands at wavelengths greater than 
9.1 « were equal for the two pellets within 2 or 3 
percent. In the range 3 to 9.1 w bands were either 
too nearly Opaque or were too close together to per- 
mit reasonable base lines to be drawn. The com- 
bination bands near 2.4 yw had transmittances of 
78.7 and 86.1 percent for the pellet containing more 
and less potassium bromide, respectively. 

In the case of alkali halide pellets the distribution 
of sample particles in the mixture is determined by 
the density of the suspending medium before press- 
ing, that is, by the volume of material plus included 
air spaces. But if m then p, appears in (7) 
only in factors m,/p, transmittance should 
practically independent of p, also. In the case of 
powder suspended in liquid the question of which 
value to use for p, does not arise. 

Figure 3 shows calculated curves plotted on semi- 
logarithmic coordinates of J,)/Z versus m, for a con- 
densed homogeneous film and for suspensions made 
with three different particle sizes. It is seen that 
for samples prepared from powders there are ap- 
parent “deviations from Beer’s law’ in the sense 
that straight lines are not obtained. The deviations 
are not great for small radii; however, it is seen again 
that apparent absorption is greatly reduced for 
larger particles or clusters. 

Because there is a finite probability that radiant 
energy will traverse a pellet without encountering a 
sample particle, it is expected that the transmission 
for opaque particles will be somewhat greater than 
the instrumental transmission zero determined with 
an opaque shutter in the optical path. This is 
equivalent to stating that J/J) will not become zero 
as a approaches infinity in eq (7). Caleulations 
with the parameters a= ©, p,»=1 gem’, 
0.714 em? m,=0.01 ¢. and m 


yam, 


sO be 


using eq (7 
p.==0.2 g/em’, S 














0.224¢ predict that for particles up to about 47 4 
in diameter the transmission zero will be rai 
than 1 percent, for 72 uw diameter particles the zero 
will be raised pereent transmission, and for 
200 uw diameter particles the zero will be raised to 34 
percent 


raised less 
to 5 


transmission 


Now the sample is to be considered as a multi 
component system. If the components are dlis- 
tributed uniformly throughout the particle, it is 
sufficient to define a in (3) by 

a S aie 


in which a; and ¢; are the absorption coefficient and 
molar concentration of the 7th component This 
reduces to a Ail 4 for Pa | frequency chosen so that 


a It is sufficient then to discuss (7 
relative to changes in a. Figure 4 contains a plot of 
Il aa, for two particle sizes and a con- 
densed homogeneous film with arbitrary values as- 
sumed for parameters m,, p,, pr, and S. Devia- 
tions from linearity are serious for very large radu 
but can probably be disregarded for small particles 
uniformly distributed in the suspension. The devia- 
tions are seen to be less serious for small a. If the 
components are not uniformly distributed through- 
out the individual particles, the calculations are more 
diffeult. An example of such a situation might 
occur when material undergoes chemical changes 
while in the form particles. In degradation 
studies, for example, degradation products might 
be found more concentrated at the surface of the 
particles than in the interior 
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Ifaisa function of p, the distance from the partiel 
center, then Beer’s law for a cylindrical increment 


Is 


of volume with radius 


dl , J ( ve 
( pid CAX\D ( ,,add ° 
] ty ‘oy | I 
The average transmittance of the sphere Is 
I] 2rl exp a(p ir od z/ar* Ih. 
p e a 


The integration is not a convenient one even for 
simple forms of a(p It will not be discussed further 
except to stress that it Is dependent on the fune- 
tional form of a It is therefore not legitimate t 
nssume a linear “Beer's law plot” for solid particles; 
undergoing chemical change 

We turn now to the case of transparent particles 
suspended in an absorbing medium. ‘This situatiol 
occurs frequently when mulls are prepared with 
mineral oil, for example It also oceurs when the; 
pellet technique is used with potassium bromide 
that contains water. The limiting behavior 
of eq (9) for vanishing particle size may be seen by 


SOTHC 


noting that [e"*(2ra—}] 1|/2e?r? approaches ¢* 
asymptotically as r approaches zero. Further sub} 
stitution of two terms of a power series for the 
exponentials leads to the expression 

lim J/J,=exp|—a(t— 3,)], (10 


in which ¢ is the thickness of the suspension and f 
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s the thickness that the powder alone would have if 
| it were in the form of a condesned homogeneous film. 
| The occurrence of the factor 3/2 in the exponent is 
) surprising. Of course as r approaches zero it Is no 
reasonable to integrate Beer’s law over a 
sphere, as was done in deriving these equations 
Values calculated from eq (9), using assumed param- 
are plotted versus particle size in figure 
| In general, samples containing the same amounts of 

material transmit more energy as the size of the 
The increase will 
occur at a greater rate for bands of higher absorption 
coefficient. The change in appearance of a hypo- 
thetical absorption band as particle size Is varied is 
demonstrated in figure 6. The envelope labelled 
“ideal” represents liquid containing no powder as 
riven by eq (10 without the factor 3/2 

While obtaining spectra of suspensions of poly- 
vinvl chloride in liquid media using a double- 
beam spectrometer, it was observed that complete 
compensation for bands arising from the liquid could 
not be obtained. For example, if a point on each 
side of the band were compensated exactly by adjust- 
ing a variable thickness cell, then that portion of the 
band falling between these points would be over- 
compensated and other portions would be under- 


? 


longer 


eters, ») 


compensated. Furthermore, when a quantity of 
compensating liquid calculated on an “ideal” basis 
1. e., from eq (10) without the factor 3/2) was 


placed in the reference beam, bands were strongly 
overcompensated. Seatter alone would not be ex- 
pected to produce these symmetric changes in band 
However, the theory involving 


shape. present 
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absorption alone gives a resonable qualitative 
description. For if an ideal-band envelope is caleu- 
lated by using a thickness appropriate for exact com- 
pensation at points on the side of the “suspension” 
band, it will have greater and less transmittance 
than the suspension band for frequencies inside and 
outside these points, respectively. The difference 
between the two contours will produce a differential 
band similar to the Attempts at 
compensation would be expected to be more success- 
ful with small particles. 


ones obsery ed. 


4. Summary 


A simple but plausible model of solid particles 
suspended in a transparent medium, with the postu- 
late that refractive index changes over molecular- 
vibration absorption bands are sufficiently weak 
to be ignored, has yielded the following predictions: 

1. The apparent absorption dependent on 
particle shape. 

2. The apparent absorption is dependent on parti- 
cle size. Samples containing large particles trans- 
mit more energy than samples with the same mass 
of small particles. Absorbance as a function of 
mass of sample for a one component system deviates 
from linearity most markedly when large particles 
are used. For smaller particles the absorbance 
curve becomes linear and in fact coincides with the 
Beer’s law plot obeyed by a condensed homogeneous 
film of the same material. In multicomponent 
systems the form of J/J) depends upon the way in 
which the components are distributed in the par- 


Is 








ticles. In the case of uniform distribution within 
the particles, deviations from linearity in absorbance 
versus concentration are greatest for large particles, 
ana agam absorbance approaches that of a condensed 
homogeneous film for small particles. All the 
marks made here regarding particle size apply to 
clustet size mM the cause of incompletely dispersed 
The desirability of obtaining small par- 
ticles and the necessity for controlling particle SIZC 
distribution and for producing complete dispersion 
of particles in the suspension medium are clear 

» The apparent absorption is not reatly affected 
by variation in and of ispending 
material, 

t. In addition, Of transparent 
particles in absorbing media has led to the conclusion 
that apparent transmittance a suspending me- 
dium should increase as the size of the party les which 
it contains 
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during the development of experimental technique | 
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Branching Ratio in the Decay of Polonium-210 
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PAMIMNA-TAy 
22 + 0.06 


' of po as found to be equal to 22 is 
ised on the assumption that the energy of the main alpha-particle group from polonium-210 

301 Mev In the calibration of the scintillation counter the angular anisotropy of the 
wma ravs from cobalt-60 was found to be 1.164 0.002. 


1. Introduction 


The number of gamma quanta emitted per dis- 
integration of polonium-210 was recently determined 
by Grace, Allen, West, and Halban [1]! in the course 
of the investigation of a method for standardizing 
polonium sources by measurement of the hard 0.804- 
\Mev [2. gamma radiation. The intensity of the 
alpha radiation emitted from a 200-me source of 
polonium-210 deposited on one side of a platinum foil 
was measured directly by means of a low-geometry 
defined-solid-angle method, whereas the O.804-Mey 
cumma-ray intensity was determined in terms of a 
source of cobalt-60 by means of a evlindrical alumi- 
num Geiger-Miller counter of unspecified filling. 
The absolute disintegration rate of the cobalt-60 
source, Which it is assumed was of a similar geometry, 
to that of polonium-210, was determined by means of 
method of beta-gamma coincidence counting 
whereas the relative efficiencies of the Geiger-Miller 
counter to the gamma radiations from cobalt-60 and 
polonium-210 were assumed to be the same as those 


determined for a counter of similar geometry by 
Bradt. Guegelot, Huber, Medicus, Preiswerk, and 
Scherrer [3]. 


The opportunity has arisen to repeat the measure- 
ment of the branching ratio of polonium-210 using, 
however, a micro-calorimetric method to determine 
the rate of alpha-particle emission and a sodium- 
odide scintillation counter, of high detection effi- 
ciency, to determine the intensity of the 0.804-Mev 
gamma rays by comparison with a cobalt-60 source, 
the gamma-ray intensity of which was determined by 
the method of gamma-gamma coincidence counting 


2. Preparation of the Polonium-210 Sources 


Experiments were performed on two polonium-210 
the first of which, however, was too weak to 
give satisfactory awecuracy This had been supplied 
sealed in a nickel tube about 1.5-em long and 0.15- 
em in external diameter. The gamma-ray spectrum 
indicated some energetic gamma-emitting contamina- 
tion, but with an activity corresponding to 3.7 me 
only a very approximate value for the branching 
ratio of 1.1 X 107° gamma quantum per alpha particle 
was obtained. 


sources, 


Figures in bracket ndicate the { tu references at the end of this paper 
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The second polonium-210 source was prepared 
from some 60 me of solution that had been care- 
fully purified prior to delivery. This purification 
consisted of depositing polonium-210 onto a silver 
foil from a stock solution containing approximately 
130 me of polonium-210. The silver foil was then 
washed free from solution and dissolved in con- 
centrated nitric acid. After expelling most of the 
nitric acid by evaporation the silver was precipi- 
tated by means of hydrochloric acid. The solution 
was filtered from the silver chloride and its acidity 
adjusted for deposition on a second silver foil, 
which was again processed in the same way as the 
first, the filtrate from the silver chloride being the 
final solution containing some 60 me of purified 
polonium-210. 

Upon receipt of this solution its acidity was ad- 
justed to a value of pH equal to 1 and 0.1 g of finely 
powdered silver was allowed to stand in the solution 
overnight. After standing for this time the solu- 
tion was filtered through a plug of glass wool in a 
glass tube having external and internal diameters 
of approximately 3 and 2 mm. The silver powder 
that was retained by the glass-wool plug was ulti- 
mately found to have extracted some 50 me of the 
polonium-210. The glass tube was then fused off 
below the glass-wool plug and above the silver 
powder to give a source 37-mm long and whose 
approximate form is illustrated in figure 1. On 
examination of its gamma-ray spectrum, this source 
was found to be free from the contamination that 
had been present in the first source, 
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FigtrRe 1 Polonium-210 source 








3. Measurement of the Alpha-Particle 
Activity 


The alpha-particle intensity of the polonium-210 
determined by means gold-cup 
radiation balance that was constructed to compare 
four national radium standards [4]. The energy 
emission was combined with an average value ? of 
5.301 Mev for the energ\ of the polonium-210 
alpha particles to give their intensity 

As this radiation balance had previously 
used primarily for relative measurements, it 
carefully recalibrated so as to give absolute rates of 
energy dissipation. Such a calibration 
essentially in determining an average of 
dk/dT for the two Peltier couples, and this 
purpose a completely new set of compensated cali- 


source was of a 


been 
Was 


CONSIStIS 
value 
for 


encased in copper sheaths and with 
Previ- 
hac 


to 


bration coils, 
permanent potential leads, was constructed 
ous calibrations had been made with coils that 
been wound on of such dimensions 
simulate the United Kingdom and the United States 
Hénigschmid primary radium standards and the 
Canadian national radium standard [4]. The re- 
sults obtained for d/id7 with the “Hénigsehmid’’, 


glass ius 


“Canadian”, and the new copper-sheathed coils 
are shown In table 1 
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The source of polonium-210, that is illustrated in 
figure 1, further important 
check to be made in connection with the operation 
of the radiation balance. It has 1] 
that the thermal “‘mixing” in the gold cups is of such 
a magnitude as to render any differences in tempera- 
ture distribution in the cups due to difference in 
source geometry unimportant. <A series of measure- 
ments was carried out therefore with the polonium- 
210 source, and a dummy source containing a glass- 
wool plug and inactive silver powder, in which the 
energy emission from the source in the erect position, 
as in figure 1, was compared with that as measured 
with the source and dummy upside down. In the 
erect position the polonium-210 is concentrated in 
the lower part of each cup in turn, whereas in the 
upside-down position a_ relatively thermal 
insulator, the glass-wool plug, is in close proximity 


Is such to enable i 
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with the Peltie: couples at the bases of the cups 
As is seen from table 2, there is no significant} 
difference between the measurements in the ty ul: 
positions; the difference between the first — tw 
measurements is not outside the possible range Ol 
experimental error | The 1 
Finally, in table 3 the results are shown for rate) from 
of energy emission from the polonium-210 > souree! ing t 
on various dates of measurement Also shown are! being 
the values that have been obtained after correctior ravs: 
using a half-life of 138.39 davs [7], to a zero date of} the o 
13.15 hours EST on June 29, 1954, on which date) obtai 
in the morning, the vamma-ray mtensitv was de 
termined 
Using an alpha-particle energy of 95.50] Ni ’ 
which corresponds to a total disintegration energ 
of 5.404 Mev, and the average value of 1276.9 4 
for the energy eCMIssion, the alpha-parti le activit 
at the above time is found to be 1,475°>< 10° disinte- 
grations per second As the gamma-ray measul 
ments were made some 3 or 4 hours earlier, it 
appropriate to correct this value to 1,476 10° dis 2 '°° 
integrations per second at the time of the gamma-ray = 
experiment 2 80 
4. Measurement of the Gamma-Ray Intensity ° ,, 
The intensity of the gamma radiation from tl 
same polonium-210 source as used in the above ” 
experiment was determined by comparison with 
cobalt-60 source, in the form of cobaltous chlorid 20 


in aqueous solution, through the use of a sodiun 
iodide seintillation high detection eff 
ereney The disintegration rate the cobalt-t 
source had been determined absolutely by the corn 
dence method ' 


counter of 


of 


— — 
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The experimental procedure was as follows: The 
polonium -210 and cobalt-60 sources of identical 
dimensions as in figure 1 were placed alternatively 
on the axis of, and 3 in. from the face of, a evlindrical 
-odium-iodide, thallium-activated, erystal, 5 in. in 
diameter and 4 in. long, as shown in figure 2. This 
erysta! was hermetically sealed in a thin aluminum 
container, using magnesium oxide a diffuse re- 
dector on all surfaces except that which was viewed 
hy a diameter photomultiplier tvpe K1198 
fhe output of the photomultiplier tube was fed 
through a linear amplifier into a single-channel differ- 
ential pulse-height analyzer from which the intensity 
versus pulse-height spectra shown in figure 3 were 
obtained As the sodium-iodide crystal is not quite 
totally absorbing, the relative efficiency for the 
detection of the gamma rays from the polonium and 
cobalt must be lated. This relative 
efficiency is given by 


is 


o-il 


SOUTCES cale 


sin BdB 


sin BdB 








(1). 
gamma ray 


efficiency, ¢, by numerical integration of eq 
The disintegration rate for the polonium 
is given by 
A (Po*") 
(Po?!) 


N ] - 1. 60 
/ , eA (Co) \ (Co), 


(2) 


where A is the area under each of the spectra of figure 
3. The area under the cobalt-60 spectrum is given 
proper weight to take into account the cascade nature 
of the gamma radiation. Four determinations of the 
relative areas give a value A (Po*°)/A (Co®°)=0.648, 
thus we obtain a value for the disintegration rate of 
the polonium gamma ray in terms of the Co®° dis- 
integration rate, 1. ¢. 
No, (Po?'® 0.5770.N, (Co). 

The estimated error of this relationship may be as 
great as 5 percent because of the uncertainty in cor- 
recting for secondary effects, such as backscattering 
into the detector. 

The disintegration rate of the cobalt-60 source 
was determined by the gamma-gamma coincidence 
method. This disintegration rate in terms of ob- 
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servable quantities is given by 


where N,; and N. are the counting rates in each of 
two similar counters, NV, is the number of coincidence 
counts between the two counters, © and © are the 
ratio of efficiencies for the two gamma rays in the 


first and second counters, respectively, and /(#) is a 
function depending on the geometry of the soure: 
and any angular correlation between the 
For the cobalt-60 source 


eounters, 


coincident gamma rays 
where there is a definite angular correlation, this 
funetion takes the form (9 
Atl, | ed: m “() 102027 Cos 0 
J J Y0009 LPs (cos 3 
ML DY |1—P, (cos 8)| / 0.1020P; (cos 3 
0.0091 [P; cos 34-2 P,(cos 3 | 
The terms containing J, are attenuation coefficients 


due to the finite angular resolution of the counters. 
where the J, is given by 


Sin pbdp » 


P, (cos By ( 


The expression for (8) ts of the same form as in eq (1 
The last term in eq (4) is a correction term [10] due 
to the finite length of the The numerical 
coefficients of the Legendre polynomials are theoret- 
ical values for a cascade transition between states with 
spins 420 [11] 

With a geometry as shown in figure 4, the theoret- 
ically expected anisotropy is— f(m)/f(2/2 L.1o0 
The experimentally observed value was 1.147, in- 


source 


dieating that when the corrections for the finite 
sizes of source and counters are accounted tor the 
measured differential anisotropy of 1.164 is well 


within the probable error (approximately | percent 
of being the same as the theoretically expected value 


of 1.167 [11] 
Several determinations of the disintegration rate 
under different reometrical conditions gave a value 


3.087 — 0.030 


~ Fe > , p : ; f ~ 
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measurements on cobalt HU source 
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Using this value, we obtain as the intensity of the | 


polonium-210 gamma rays 
j 
Vo, (Po 


1.782 5Y, lO? dps 


— 


5. Branching Ratio | 


Combining the gamma-ray and alpha-partiel 
intensities, the following value for their ratio 
obtained 

Voy _1.782¢ 10! 

Vow 1.47610 

[ 

1 .20- 0.06 107’ gamma quantum pet alpl 
particle 

This value of (1.21 + 0.06) 107° is to be compared | 
with that of (1.8+0.14 10 obtained by Grae 
Allen, West, and Halban [1] and (1.6402 10 
by M. Riou [12]. The internal conversion 9 
K-shell electrons in’ polonium-210 is 0.0084 
and correction for L-shell conversion electrons 14 
vives a total conversion coefficient of 0.0107 Using, 
this value, the branching ratio is found to be 1.22 


' 
| 


10 Dd percent 
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Mechanism of High-Speed-Waterdrop Erosion of 
Methyl Methacrylate Plastic 


Olive G. Engel 


\ mechanism for the high-speed-waterdrop erosion of methyl methacrylate plastic is 
postulated as a result of a parallel study of the damage that is done by the impacts of steel 
-pheres and of deforming lead pellets. tesults of the parallel study are described. The 
mechanism can be extrapolated with safety only to other brittle materials of low tensile 
strength, the other characteristic properties of which are closely similar to those of methv] 


methacrviate plastic 


1. Introduction 


that 
surfaces at 


The waterdrops cause when they 
strike solid very high velocities has 
recently received considerable attention because of 
the damage that is produced when high-speed air- 
craft fly through rain. This erosion is essentially the 
same as that which occurs in steam turbines because 
of the waterdrops in wet steam. The process has 
been studied for the past 25 vears [1],! but 
understanding of the complete damage mechanism 


erosion 


no clear 


has emerged 

{nv attempt to explain the damage that results 
Irom high-speed-waterdrop impingement must take 
into account the properties of the waterdrop under 
impact conditions because the high-speed waterdrop 
is the damage tool. That is, the destructive force 
that causes this type of damage is a result of the 
impact of the waterdrop against the solid surface. 
Under impact at a velocity of the order of the speed 
of sound in air, water acts as though it were hard, but 
it retains its liquid property of flow. It has recently 
estimated [2] that the Impact that 
results from a collision of this tvpe at an impact 
velocity of 600 mph is of the order of 30,000. psi 
and that the radial-flow 1,400 
mph 

The damage that results from these high-speed 
liquid-to-solid collisions is equally a function of the 
properties of the structural material that sustains 
the damage The response of materials of different 
properties will be different enough to introduce 
notable modifications in the erosion process. 

The analysis of a difficult problem may sometimes 
be accomplished by reducing it to a simpler problem 
that is easier to understand. It was suggested by 
Schréter 3] in recard to cavitation for 
example, that physically parallel studies might’ be 
rewarding than further observations of the 
cavitation progress itself because the cavitation can- 
not be reduced to slow motion nor separated from 
other phenomena. As just pointed out, under im- 
pact conditions a waterdrop behaves like a hard 
sphere, but unlike a sphere of hard, solid material, it 
retains its liquid property of flow. The results of 
the parallel! study of the Impact of steel spheres, and 
of lead pellets, on methyl 
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erosion, 


more 
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methacrylate plastic appears to indicate the mech- 
anism of high-speed-rain erosion. 


2. Collision of Steel Spheres and of Deform- 
ing Lead Pellets with Methyl Methacrylate 
Plastic 


The plastic blocks used for this study were of 
Plexiglas I-A obtained from the Rohm & Haas Co. 
in \4-, l-,and 3-in. thicknesses. The projectiles were 
32-In. steel spheres, 0.22-in.-diameter soft deforming 
lead pellets that were shot with small charges of 
powder from a 0.22-caliber gun, and %¢-in. steel 
spheres that were shot from a 0.30-caliber gun. The 
impact velocities were determined by momentum 
transfer to a pendulum through the block of plastic. 

The impacts resulted in a characteristic type of 
fracture of the plastic. At relatively low velocities 
of impact, cireles of short, very fine, straight cracks 
were observed (see fies. 1, 2. and 3). These straight 
cracks were perpendicular to radii from the central 
point of impact. The center spot of the impact was, 
however, completely without any such cracks. At 
relatively high velocities of impact, the circular 
region of fine cracking on the impact surface was 
bounded by avery well-defined circular crack when 
the projectiles were steel spheres. Such a crack was 
not observed at these velocities when the projectiles 
were lead pellets. However, in the latter case, a 
circular subsurface shear crack was produced. These 
subsurface shear cracks appear as large dark circles, 
with their center at the center of the damage site in 
figure 3. They are out of focus with the surface 
details. 

Radial cracking resulted from the impacts of 
steel spheres but not from those of deforming lead 
pellets. On thin plates of Plexiglas I-A and at low 
velocities the radial cracking appeared first on the 
side of the plate opposite to the impact face. At 
high velocities on these thin plates the radial cracks 
penetrated through the plate to the impact surface. 
On the thick plates, very short radial cracks appeared 
on the impact surface itself. They were observed 
only at relatively high velocities. 

A detailed discussion of the specific types of crack- 
ing and of a model that explains the observations is 
contained in the following sections. 
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PIGURI Damage site produced on Plexiglas I-A as a result 
amu lip e shot with de forming lead pe llets. 
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2.1. Circular Crack and Fine Cracking in the 
Circular Zone 


At low impact velocities both for impacts of steel 
spheres and for impacts of deforming lead pellets the 
damage consisted only of a few circles of fine, short, 
straight cracks that enclosed the central point of the 
impact. The central point of the impact itself was 
free of anv cracking 1, 2, and 3 As the 

| velocity was increased, the amount of fine cracking 
| increased, and the diameter of the zone of fine crack 
| ing increased. Finally, a velocity was reached at 
} which a well-defined crack formed around the cir- 
' cular zone of fine cracking when steel spheres were 
This outer crack was 
too deep for a depth measurement with the inter- 
lerometer. When deforming lead pellets were used, 
subsurface shear cracks formed at high velocities. 
The subsurface shear cracks appear as broadened 
ares Of a circle out of focus with the surface details 


see figs. 


used as projectiles see fic’, | 


see flo’. 3). 

The diameter of the circular zone of fine cracking 
lor comparable velocities was influenced by the 
thickness of the plastic plate only in the case where 
the plate thickness was as low as 4 in. Comparison 
of the diameters of damage sites on 1-in.-thick and 
on 3-in.-thick of Plexiglas I-A for %¢-in. 
| spheres at comparable velocities showed that they 


blocks 


53 


FIGURE 4 C 
-defined cl 


O88 ection of a 


region of fine cracking bounded 


hy a well cular crack similar to that shown in 


figure ] 


were essentially the same. Comparison of the diam- 
eters of damage sites on \-in.-thick and on 3-in.- 
thick Plexiglas I-A for %s-in. spheres at comparable 
velocities showed that the diameters on the \-in.- 
thick material were larger. 

The diameter of the circular zone of fine cracking 
for comparable velocities also varied with the diam- 
eter of the steel sphere. Comparison of the diameter 
of the damage site produced by a %-in. sphere with 
that produced by a %¢-in. sphere on \-in.-thick 
Plexiglas I-A at the same velocity showed that the 
damage site produced by the larger sphere was the 
larger. 

An interferometric study of the circular region 
bounded by fine cracking indicated that the cross 
section of this region must be as indicated in figure 4. 

When deforming lead pellets were used, the cracks 
were widened and material was broken out of the 
surface. The widening of the cracks was in the 
direction of the radial flow of the lead. The break- 
ing out of material was along cracks, was in the 
direction of the flow of the lead, and was on the side 
of the crack away from the center of the impact (see 
figs. 2 and 3). The widening of cracks and the 
breaking out of material was notably absent where 
steel spheres were used as projectiles. 


2.2. Radial Cracking 


Cracks appeared that went out radially from the 
circular crack in the case of the impacts of steel 
spheres. They were especially evident in the dam- 
age marks left by the impact of steel spheres against 
the thick blocks of plastic. These cracks extended 
only a very short distance within the circular region 
of fine cracking. In figure 1 they appear as feather- 
like structures around the circular crack. 

There was also another type of cracking that 
occurred as a result of the impact of steel spheres 
when the plastic plates were thin. This was a 
radial, or star-shaped, cracking similar to that which 
results from the impact of a small, hard object 
against thin glass plates. In the case of \-in.-thick 
Plexiglas I-A these radial cracks appeared on the 
side of the plastic opposite the impact face at rela- 
tively low velocities. They converged under the 
central point of the collision. They did not extend 
all the way through the plate to the impact surface. 
Furthermore, there was no circular crack or circular 
area of fine cracking on this reverse face of the plate. 
At high velocities these radial cracks extended com- 
pletely through the plate to the impact surface. 

These radial cracks were quite long. In the case of 
impacts of %.-in. spheres at the relatively low 
velocities of 232 and 269 ft/sec the cracks were over 
2 in. long and for the higher of these two velocities 








were close to 3 in. lone. On the other hand, radial 
cracks that resulted from the impact of 7o-in. spheres 
at the high velocities of 611 and 943 ft/sec were 
only about 1% in. long. From this it would appear 
that the length of these cracks may be a function 
both of velocity and of sphere diameter. It would 
that the total crack length should correlate 
with the energy transferred to the plastic by the 
sphere 

Bowman, Smith, and Kies [4] have plotted total 
length of radial cracking against velocity and have 
found a maximum in the curve for Plexiglas I-A 
at a velocity of 775 ft/sec, and for Plexiglas LI] at 
a velocity of about 875 ft/see They were using 
*e-In.-diam steel spheres The shots were made 
against 6-in.-square plates of plastic that were '; in 
thick. The reported maximum is a suppressed maxi- 
mum that was observed under these experimental 
conditions. They found, furthermore, that Plexiglas 
I] that had been subjected to hot-work reached 


seem 


maximum crack length at a lower velocity than 
Plexiglas II in the as-received condition 
2.3. Crazing in a Circle 

The fine cracking that has been observed at the 

damage sites produced by the impacts of steel 


spheres and of deforming lead pellets on Plexiglas 
I-A may be designated as crazing. In the case of 
the cup-shaped depression that must form when 
a steel sphere or lead pellet strikes a plastic surface, 
the tensile stresses extend as radii from the central 
point of impact, which is under compression, to the 
rim of the cup-shaped cavity. The fine cracking 
that forms a result of these tensile i 
more or less parallel in tiers that form at an angle 
another, but the individual 
not mutually paralle! 

Recently, Wolock and Newman [5] traced 
the cross section of a craze crack of the generally 
accepted type with the interferometer. They have 


nS stresses 1s 


to ome fine cracks ure 


have 


found that the surface material is raised on both 
sides of the crack (see fig. 5 It may be postulated 
that such cracks form as stress relievers after the 


material has reached the limit of stress relief from 
plastic flow, or if plastic flow cannot take place fast 
enough. If plastic flow has occurred either prior 
concomitant with the formation of the craze 
the material will tend to come back into its 
original position after the stress is relieved, and the 


to o! 


crack, 


crack will be raised on each side as is observed. As 
a part of the present investigation, Dr. Newman 
made a similar interferometric study of the fine 
cracking that was produced by the impact of a steel 
sphere on Plexiglas | A and found that a raised 
rida exists on both sides of the crack. as shown 
in fieure 5 
Y“ i 
ad 

| ' ) ( a / f j ‘ 

nad he nile 


2.4. A Model for the Impact Damage of Deformin, 
Lead Pellets and Steel Spheres on Methyl Methg, 
crylate Plastic 


The observations that have been described eq 
be explained rather easily by a simple model } 
which it is assumed that the impact of a steel spher 
or of a deforming lead pellet produces a cup-shape 
cavity in the plastic by compressing it.  Tengil 
should then be introduced in’ the 
indicated as A, B, and C in figure 6. The stresses 
in region A should become important as soon as any 
depression takes place, regardless of how: shalloy 
if may On the other hand, that ar 
introduced in the regions indicated as B and C may 
only become appreciable after the cup-shaped cavity 


stresses regions 


stresses 


be 


| 


| 


~~ 


—— 


reaches the more acute stage of depression indicated| 


schematically in figure 6. The stresses in region | 
exist over the relatively sharp angle between th 
undisturbed surface and the wall of the depression 
It is postulated in section 2.3 that craze cracks 
may form as stress relievers after a material that 


under tensile stress has reached the limit of 


‘ 
/ 


Stress 


relief from plastic flow or if plastic flow cannot take! 


place fast enough. Similarly, the fine cracking 


which appears to be crazing, should serve to relieve) 


the tensile stresses, indicated as A in figure 6, 
are produced as the plastic material is depressed by 
the impact A cavity as deep as that indicate 
schematically in figure 6 will only be produced at 
high impact velocities and only by a sphere thai 
does not flow on impact. Hence, the tensile stresse 
at B should only become acute at high velocities 
and then only in the case where a steel sphere rath 
than a deforming lead pellet is the projectile. Thes 
may be responsible for the well-define 


stresses 


thar' 


circular crack that forms at high impact velocitie| 


around the circular region of fine cracking in tl 
damage sites produced by the impact of steel spheres 
The tensile stresses, indicated as C in figure 6, may 
be responsible for the feather-like eracks that vO Ol 
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radially from the circular crack in the case of the 
impacts of steel spheres. 

On this model, figure 6, there should be severe 
compression during the impact where the hard 
sphere 1s bearing against the plastic to form the cup- 
shaped cavity. MeAdam, Geil, and Jenkins [6] 
have pointed out that material under compression 
does not tend to undergo a tensile failure. It 
notable that there is no fine cracking in the severely 
compressed region at the center of the damage sites 
see figs. 1, 2, and 3 

In the light of the model that has been assumed, 
essentiaily the same diameter of circular cracking 
would be expected to occur on a 1-in.-thick plate as 
on a 3-in.-thick plate of plastic, that is, for those 
plate thicknesses that are large in comparison with 
the shot diameter. It is also expected that for a 
much thinner plate, the diameter of the damage site 
would be increased by the bowing of the plate as a 
whole around the impinging sphere or pellet. This 
was observed to be the case 


Is 


Furthermore, on the model that has been assumed, 
a gradient of tensile stress would be expected up the 
sides of the cup-shaped depression, that is, along the 
arrow marked A in firure 6. Griffith \7] derived an 
expression that showed that there should u 
maximum crack length for every value of stress that 
produced a crack. A magnification of the fine- 
eracked circular region produced by one of the shots 
shows some evidence that the lengths of the fine 
eracks at any given distance from the central point 
of impact are, in general, qualitatively of the same 
size, and that there is a size gradient from the rim 
to the bottom of the depression. However, there 
are mans exceptions to this observation. 

The radial cracks that are observed when the 
projectiles are steel spheres and that appear first on 
the reverse side of the plate ean be explained as 
tensile failures. After the compressional wave pro- 
duced by the impact on a relatively thin plate (one- 
dimensional case) has just moved through the plastic, 
the plastic material on the free face opposite the 
impact surface should tend to bulge when the com- 
pressional wave reflects there as a tensile wave. It 
is possible that a star-shaped hole may open to 
relieve the After the compressional 
wave has reflected as a tensile wave, the plastic 
material should snap back into a star consisting of 
only A search was made for this 
effect by Hl. L. Smith [8], using a high-speed camera 
operating at 15,000 frames a The star- 
shaped crack appeared, however, in the space of one 
fram It is possible that the stress may be relieved 
sufficiently by crack formation alone, and that a 
star-shaped hole may never form. A camera having 
a higher operating speed will be needed to determine 
this 

The deforming 
flowed on impact 


be 


tensile stress 


radia! cracks 


second 


lead pellets used as projectiles 
to approximately twice their 
original diameter. This flow of the lead produces a 
loree against any surface irregularity in its path 
which, in conjunction with the restraint of the ir- 
regularity by the underlying material, results in a 


55 


torque tending to break the material in the irregu- 
larity away from the underlying material. The 
raised ridges along the fine cracks could constitute 
bearing points for this torque and offer an explana- 
tion for the observed fact that where lead pellets 
were used as projectiles the cracks were widened and 
material was broken out of the surface along the 
cracks and in a direction away from the center of 
the impact (the direction of the flow of the lead) 
2 and 3). 

This flow property of the lead also explains the 
subsurface shear cracks that appear as large dark 
circles out of focus with the surface details in figure 
3. The flow of the lead over the surface while it is 
simultaneously compressed against the surface would 
tend to move the surface layer of the plastic with 
respect to the under lavers. The subsurface cracks 
are not observed in the damage sites produced by 
the steel spheres, which did not flow as a result of the 
collision. 

It is the flow property of the deforming lead pellets 
that makes the damage sites produced by them dif- 
ferent from the damage sites produced by the steel 
spheres. The fracture that results from the impact 
of deforming lead pellets is more nearly related to 
what must happen when liquid drops strike a sur- 
face than is the fracture that results from the impact 
of steel spheres. In the case of waterdrops, of 
course, the damage from a single blow cannot be as 
decisive as that which results from a single impact 
of a deforming lead pellet because the density of 
water | much lower than that of lead. The 
damage caused by high-speed impact of waterdrops 
on methyl methacrylate plastic is discussed in detail 
in the following sections. 


(see figs. 


Is so 


3. High-Speed-Rain Erosion of Methyl Meth- 


acrylate Plastic 


The preceding study of the impact of steel spheres 
and of deforming lead pellets on methyl methacrylate 
plastic is a physically parallel study to the high- 
speed-rain erosion of brittle plastics of low tensile 
strength. To make a realistic approach to the 
mechanism of the high-speed-rain erosion process, it 
is essential to keep in mind that the projectile is a 
water sphere that acts as though it were hard but 
that is capable of an extremely fast radial flow the 
instant after it impinges. 

The process of high-speed-rain erosion should be 
very similar to the process of producing damage by 
impacts with deforming lead pellets. For plastics 
that are similar in properties to methyl methacrylate 
plastic, radial fine cracking should form around the 
point of impact of a hard, impinging, spherical water- 
drop due to the tensile stresses that are produced 
when the compressive impact stress forms a cup- 
shaped cavity in the material. A breaking out of 
material from the surface should be observed as a 
result of the rapid radial flow of the water after the 
collision over the raised ridges of the fine cracks that 
have formed 
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3.1, Circular Fine Cracking as a Result of High- 
Speed-Water Blows 


A search for damage sites comparable to those that 
are observed as a result of impacts of deforming lead 
pellets was at first unrewarded. Specimens of 
methyl methacrylate plastic that were rotated at a 
velocity of 500 mph through 1-in./hr artificial rain 
at the Cornell Aeronautical Laboratory for periods 
of 30 to 40 sec were already in an advanced stage of 
erosion. The first damage sites were obscured by 
later overlapping impacts, so that nothing could be 
concluded in regard to their initial configuration. 
No gene ral pattern was \ isible, and all that could be 
concluded from them was that the erosion damage 


went out from cracks and crack intersections that 
may or may not have been there before the water 
struck. 


However, individual damage sites could be seen on 
aspecimen similarly exposed for only 25 see. Figure 
7 shows the result. The configurations that stand 
out in the photograph are the cireles of fine cracking. 
Here again, the circles of parallel lines in tiers can be 
seen. It is also notable that here again the craze 
cracks meet at an angle. In figure 7 the widening of 


these cracks due to the outward flow of the water 
from the central point of impact can be seen. Direct 
comparison should be made with figures 2 and 3, 


where the same configuration of lines and the same 
widening of cracks can be seen as the result of the 
impact and flow of deforming lead pellets 

Figure 8 higher magnification the most 
prominent circle shown in figure 7. In figure 8 the 
angular meeting of the parallel lines is apparent as is 
also the widening of the cracks by breaking out of the 
surface material along them in a direction away from 
the central point of impact as a result of the radial 
flow of the water. t,o and & the 


of 


is a 


Krom figures 2, 3, 7, 

conclusion ean be drawn that the mechanisms of dam- 

(life wroduced hy the im wictl of de forming lead vellets 
| by pact mn } 

and by the high-speed impinge ment of waterdrops on 
ethyl methacrylate plastic aire ide ntieal. The differ- 
ce between them is one of degree only. 


3.2. Comparison of the Destructive Action of the 
Impact of Lead Pellets and the Impact of Water- 
drops 


The destructive action of these flowing projectiles 
sdue (A) to the impact of the hard sphere that con- 
stitutes the projectile itself, and (B) the radial 
flow of the projectile after the collision, as follows: 

A) The first of these two tvpes of damage is the 
production of a circle fine cracks With respect 
to this type of damage, the impacts of lead pellets 
ire more destructive than the impacts of waterdrops, 
by comparing the amount of fine 
cracks produced by a single impact in figures 2 and 3 
with that produced by a single impact in figure 7. 
The impact velocity of the lead pellets that produced 


to 


seen 


as can be 


| the damage sites shown in 


) 


figure 3 was about half 
that of the waterdrops that ‘produced the damage 
sites shown in figure 7. The lead-pellet velocity that 
produced the damage sites shown in figure 3 was 
appreciably higher than that which produced the 
damage site shown in figure 2. It is easy to see that 
the extent of damage observed in figures 2 and 7 is a 
logical result of the kind of flowing projectile that 
produced it. The impact pressure that results when 
a spherical liquid projectile strikes a flat, solid 
surface that has an infinite modulus of elasticity is 


given by [2] 


P=(a/2)cpvo (1) 
where P is the impact pressure, a is a coefficient less 
than unity that arises from the flow property of the 
liquid projectile, ¢ is the speed of sound in the liquid 
projectile, p is the density of the liquid projectile, 


and vy is the impact velocity. Lead has a much 
higher density than water and a value of @ that is 
closer to unity. Therefore, according to eq (1), a 


lead pellet should exert a greater impact stress than 
2 water sphere. The lead pellet should form a 
deeper depression than the waterdrop does in the 
surface of the plastic result of its collision, 
should produce more tensile stress in the surface 
about the depression that it produces, and conse- 
quently should produce more fine cracking. 

(B)} The second type of damage is (a) the widening 
of the surface cracks and the breaking out of material 
from the surface, both of which are observed in the 
damage sites produced by waterdrops as well as in 


as a 


those produced by lead pellets, and (b) the sub- 
surface shear cracks that are observed only in the 
damage sites produced by the lead pellets. The 


damage (a) may result from the application of the 
external torque exerted by the rapid radial flow along 
the surface against irregularities in the surface. In 
this connection the raised ridges of the craze, or fine 


cracks, may be a factor. The damage (b) appears to 
be due to a shear stress that tends to move the 
surface layer of the plastic solid with respect to the 


underlying lavers. 

Whether the damage caused by the flow as a result 
of impact of a waterdrop or of a lead pellet is con- 
sidered more severe depends on the criterion that is 
taken for the extent of damage. The waterdrop flows 


out to many times its original diameter and there- 
fore breaks a large amount of material from the sur- 
face. The lead pellet flows only to about twice its 


original diameter and hence is less effective than water 
in this respect. On the other hand, the shea 
exerted by the flow of the lead is much greater than 
that exerted by the flow of the water, because it pro- 
duces subsurface shear cracks, which the water does 
not produce even when its impact velocity is about 
double that of the lead. On the basis of its shear 
therefore, the radial flow of the lead is more 
destructive than that of the water. 


r stress 


stress, 
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3.3. Shear Stress Due to the Radial Flow of an 
Impinging Waterdrop 


The shear stress, 7, between layers of liquid in 
laminar flow is given by the product of the viscosity, 
uw, and the velocity gradient through the moving 
sheet of liquid perpendicular to its direction of flow. 
That is, 


or 
Oz 


kh 


where v is the velocity of the sheet of liquid, and Is 
the direction through the thickness of the liquid 
sheet The laver of liquid molecules in direct con- 
tact with the solid has zero velocity, but the velocity 
eradient is not zero, and the shear stress is applied 
to the solid 

Faust [9] and Hyde [10] have found that the vis- 
cosity of liquids not only increases with pressure, 
but it increases at an increasing rate as the pressure 
Hence, it seems entirely reasonable to con- 
clude that water impacts will produce subsurface 
shear cracks at sufficiently high impact velocities, 
even though thev are not observed in the views 
shown in figures 7 and 8, where the relative impact 
velocity was 500 mph. 

[t is possible to estimate roughly the shear stress 
exerted by the radial flow of waterdrop that 
strikes a solid surface at a relative velocity of 500 
mph. In developing an equation for the impact 
pressure under an impinging water sphere {2}, it 
was shown that the thickness, 6, of the sheet of 
water in radial flow when the pressure is at a maxi- 


rises, 


al 


mum is given by 


6 47( 1—a)* ( ; ). | 


is the radius of the flattened water sphere, 
a a coefficient that arises from the flow property 
of the liquid and that is very close to unity at high 
impact is the relative Impact velocity 
between the water sphere and the surface that it 


where 


Is 


\ clocities. 


strikes, and ¢ is the speed of sound in water. The 
median radius of 1-in./hr raindrops, such as were 
used in producing the damage sites in figures 7 


and 8, is 0.095 em, so that the radius of the flattened 


sphere is 9.5 em If the coefficient @ is taken to 
be 0.99, then 6 is found to be 0.00008 em 
The average velocity of the radial water flow at 
the time of maximum pressure also was found [2] to 
be eve a by 
V CA». | 


Kor 


a relative impact velocity, of 500 mph, the 

velocity, is found to be 58,000 
\s this an average value of radial 
flow velocity, a rough estimate of the velocity eradi- 
ent through the radially moving sheet of water ean 
be made by taking it to be equal to Av/Az, where A 
is half the value of 6. To this degree of approxima- 
tion, ris found to be 210 psi, | 


flow about 


Is the 


eCrhiysec 


because the viscosity 
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of water at room temperature is about 0.01 pois 
This value is too low because of the increase of 
viscosity with pressure. From the data of Faust 
the viscosity of ethyl aleohol increases by a factor ; 
about three for an increase of pressure comparab]\: 
to the impact pressure that is produced when 
waterdrop strikes a solid surface at a relative 
velocity of 500 mph. If the increase in the viscosity 
of water, which is also an associated liquid, is about 
the same as that of ethyl aleohol, then 7 is about 
630 psi for the radial flow of a waterdrop that strikes 
a solid surface at the relative velocity of 500 mph 
This is below the static shear strength of methyl] 
methacrylate plastic, which is about 8,000 psi [11]. 

The ealeulated value of 7 is rather sensitive to the 
value that The value of 0.99 for a 
was used in place of the value 0.9, which has been as- 
sumed to be the value of ain previous calculations for 
this magnitude of the relative impact velocity in order 
to make theshearing stress as high as possible. It will 
be interesting to know at what impact velocity sub- 
surface shear cracks will be found to be produced (X- 
perimentally on methyl methaervlate plastic by water 
blows. This may come to light when a method to 
produce higher relative impact velocities between a 
water sphere and a solid surface 1s develope 


is chosen for a. 


4. Mechanism of High-Speed-Rain Erosion 
on Methyl Methacrylate Plastic 


The effeet of high-speed-waterdrop impingement 
on methyl methacrylate plastic correlates with the 
damage produced in the parallel study of impacts 
of steel spheres and especially of the impacts Ol 
deforming lead pellets. In the parallel study it was 
assumed that a cup-shaped cavity is produced in 
the plastic surface as a result of the compressiv 
load exerted by the impinging projectile. It was 
further assumed that depressing the surface in this 


wav introduces tensile stresses that, for a_ brittle 
material of low tensile strength, produce cracks 
In the ease of a flowing projectile the cracks are 


widened, and material is broken out of the surface 
along the cracks and at crack intersections by the 
torque that the ultrarapid radial flow exerts 

At high impact velocities a waterdrop bel 
a hard sphere, but unlike a sphere of hard solid 
material, it retains its liquid) property of flow 
These properties of the waterdrop unde 
conditions make it comparable to the steel spher 
and to the deforming lead pellet It is logical that 
the mechanism by which it damaged the plasti 
surface should be entirely comparable to that by 
which the steel spheres and the deforming lead pellets 
damaged the surface 

Beal and Wahl [12] have found that for a 
relative impact velocity the impingement of water- 
drops at a 3-in./hr rain rate, in which the medial 
2.5-mm diam, produces more damag 


i 
aves LK 


Impact 


oivel 


drop size is 2 
than the impingment of waterdrops at a 1-in. hr rail 
rate, in which the median drop size is 1.9-iim cian 
They also found that where three nozzles, each ol 


which produced 1-inwhr rain with a median droj 
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we of 1.9-mm diam, were used, the 3-in./hr volume 
of main Was not proportionately more severe than the 
1 ,./nbr volume of rain if the drop size is the same in 

ch. They concluded that the greater damage 
»-oduced by the 3-in./hr rain rate with the median 
drop size of 2.5-mm diam must be due to the larger 
drop size rather than to the increased number of 
drops per unit time. In the light of the mechanism 
that has been postulated, it would follow that this 
should be the case Comparison of the diameters of 
damage sites produced by a %e-in. and by a %¢-in. 
steel sphere at the same relative impact velocity and 
mn methyl methacrylate plastic of the same thickness 
showed that the damage site produced by the larger 
sphere was the larger (see section 2.1). Similarly, 
the damage site produced by the average 3-in./hr 
raindrop will be larger than the damage site produced 
by the average 1-in./hr raindrop. The 3-in./hr 
raindrop will also provide a larger volume of water to 
produce damage in the process of radial flow. It is 
notable that eq (1) for the pressure produced by the 
impact of a water sphere against a flat surface is 
not a function of the drop size. 

Although it would appear that the rain-erosion 
mechanism on methyl methacrylate plastic has been 
explained by this study, one cannot safely extrapo- 
late the result to a// other brittle plastics of low 
tensile strength. As pointed out in section 1, the 
erosion mechanism is not only a function of the 
impact properties of the waterdrop, but it is also a 
function of the specific properties of the solid material 
against which the waterdrop impinges. The mecha- 
nism that has been postulated for methyl metha- 
eryvlate plastic can be extrapolated with safety only 
to other brittle plastics of low tensile strength having 
other characteristic properties that are comparable 
to those of this plastic. Furth and Born [13] have 
advanced the concept that breaking is a form of 
nelting: “melting being nothing else than a breaking 
due to the action of the heat movement of the atoms; 
" putting it the other way round, breaking is nothing 
than melting enforced by the of the 


else action 
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(see also references [14, 15, 16]). If this 
concept is correct, a brittle plastic of low tensile 
strength that differs appreciably from methyl 
methacrylate plastic either in heat conductivity or 
in actual melting temperature could respond in a 
very different manner to the high-speed-water blows. 
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